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ABSTRACT 



Nikodem J. Poplawski 



A Relativistic Description of Hadronic Decays of the Exotic Meson 7Ti 



Exotic mesons are striking predictions of quantum chromodynamics that go beyond the 
quark model. They can provide great insight into understanding phenomena such as 
asymptotic freedom, confinement, and dynamical symmetry breaking. This work analyzes 
hadronic decays of exotic mesons, with a focus on the lightest one, the J PC = 1~ + 7Ti, in a 
fully relativistic formalism. The relativistic spin wave functions of normal and exotic mesons 
are constructed based on unitary representations of the Poincare group. The radial wave 
functions are obtained from phenomenological considerations of the mass operator. We find 
that fully relativistic results using Wigner rotations differ significantly from nonrelativistic 
ones. Moreover, the S + P selection rule is also satisfied in relativistic formalism. Final 
state interactions do not change these results much. 
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Chapter 1 



Introduction 



Atomic nuclei consist of nucleons (protons and neutrons) which are built of quarks and 
gluons. Quarks and gluons can also combine to form matter called mesons that are not 
commonly found on Earth, but are naturally created in processes that occur in outer space. 
Quarks and gluons interact with each other via the strong nuclear force to form hadrons 
such as nucleons and mesons. This strong force confines quarks and gluons inside hadrons. 
Exchange of light mesons, in particular it, p, a, u, can be used to approximate the effective 
nucleon-nucleon interaction that binds protons and neutrons in atomic nuclei. In the so- 
called constituent quark model which describes matter to a good approximation, mesons 
are regarded as quark-antiquark pairs. In the same model, the proton, the neutron, and 
other baryons are triplets of quarks. 

Recently, there has been evidence for a new kind of particles dubbed "exotic mesons". 
They are called exotic because they have unusual quantum numbers, which are not allowed 
for quark-antiquark pairs. A theoretical description of the mechanism of exotic meson decays 
is vital to our understanding of nature of the quark-gluon interaction. A complete model 
which describes the behavior of exotic mesons should be based on the theory of relativity, 
in which one deals with velocities close to speed of light. In this work we will attempt to 
construct a relativistic model of hadronic decays of exotic mesons, in particular the so-called 
it i exotic meson. The main goal is to determine how much relativistic description of the n\ 
decays differ from the existing nonrelativistic predictions. 

This thesis is organized as follows. A brief review of QCD and the constituent quark 



model will be given in the following sections. In Chapter 2 we will introduce the exotic 
mesons and review the experimental situation. In Chapter 3 we will discuss the foundations 
of our model. In Chapter 4 we will present relativistic dynamics for noninteracting and 
interacting particles, and prepare the spinor framework for relativistic exotic meson decays. 
In Chapter 5 and 6, a relativistic construction of normal and exotic meson wave function 
will be given, with a particular emphasis on relativistic effects, including spin-orbit coupling 
and phase space modification. Decays of normal mesons and of the 7i"i exotic meson will be 
analyzed in Chapter 7 and 8, respectively. In Chapter 9 we will explore the effects caused 
by residual interactions between the decay products, the so-called final state interactions. 
The final conclusions will be summarized in Chapter 10. 

1.1 QCD and the constituent quark model 

The discovery of the pion in 1947 helped to understand the nature of the nucleon-nucleon 
force. However many other mesons and baryons were found shortly after, which implied 
that none of these particles were elementary, and that pions were not the quanta of the 
strong interaction. In order to extend the N + it scheme to other hadrons and to account 
for certain decay patterns such as a long lifetime of the £~, M. Gell-Mann, T. Nakano and 
K. Nishijima independently proposed in 1953 the concept of strangeness. Strong decays with 
a short lifetime on the order of 10 -24 — 10~ 22 s had the property of conserving strangeness, 
whereas much longer weak decays violated this conservation. 

In 1961, Gell-Mann and independently Y. Ne'eman introduced the eightfold way, i.e., the 
SU(3) symmetry ordering of all subatomic particles analogous to the ordering of the chemical 
elements in the periodic table. It was a generalization of the SU(2) isospin symmetry of the 
nucleon into an SU(3) with strangeness as a second additive quantum number. A success 
of this theory was the discovery of the Q~ baryon. The lightest mesons were also organized 
into a nonet. 

As in the periodic table, a large number of hadrons suggested the existence of sub- 
structure. SU(3) symmetry and, in particular, its breaking led Gell-Mann and G. Zweig to 
postulate the quark in 1963 P3E]. They suggested that mesons and baryons are composites 



of quarks or antiquarks having three flavors: u, d and the heavier s. Since fractional charges 
have never been observed, the introduction of quarks was treated more as a mathematical 
explanation of flavor patterns than as a postulate of an actual physical model. In 1965 
O. W. Greenberg, M. Y. Han and Y. Nambu introduced the quark property of color charge 
in order to remedy a statistics problem in constructing the A ++ wave function. All observed 
hadrons had to be neutral singlets of the color SU(3) symmetry. 

In 1968-69, an experiment at SLAC in which electrons were scattered off protons (deep 
inelastic scattering) led J. Bjorken and R. P. Feynman to realize that the data could be 
explained as evidence of small hard cores inside the proton called partons. This picture, 
however, had several problems, for example that ~ 50% of the proton momentum was not 
in quarks. 

In 1973, a quantum field theory of the strong interaction was formulated by H. Fritzsch 
and Gell-Mann, based on the Yang-Mills color SU(3) nonabelian gauge symmetry which 
is different from the approximate flavor SU(3). In this theory, called quantum chromody- 
namics (QCD), quarks and massless gluons (quanta of the strong-interaction field) carry a 
color charge. The structure of QCD is similar to quantum electrodynamics, based on the 
U(l) symmetry, but much richer. Because gluons carry color charge they can interact with 
other gluons. The gauge-invariant Lagrangian, unlike that of QED, has cubic and quartic 
terms in the field potential leading to nonlinear classical equations of motion and interesting 
topological properties of the vacuum. 

In 1973, D. Politzer, D. Gross and F. Wilczek discovered that QCD has a special property 
called asymptotic freedom, i.e., at short distances the coupling constant is small enough for 
perturbation theory to be valid. Unfortunately, at larger distances the coupling constant 
is on the order of 1 and the QCD Hamiltonian cannot be solved perturbatively. Quantum 
chromodynamics exhibits at this scale another distinct feature, quark color confinement, so 
that we may observe only colorless particles. At present we know six quark flavors: u, d, s, 
c, b and t. Together with gluons they are included in the Standard Model of fundamental 
particles and interactions. 

The QCD Lagrangian is in practice very difficult to solve. Thus, one is forced to deal 
with phenomenological models. One such model, the QCD sum-rule approach, was intro- 



duced in the late 1970's and applied to describe mesonic properties [3]. This technique was 
also extended to baryons |1J. The basic idea of QCD sum rules is to match a QCD descrip- 
tion of an appropriate momentum-space correlation function with a phenomenological one, 
and establish a correspondence between hadronic and quark degrees of freedom 5 . This 
approach provides a connection between the QCD Lagrangian and hadron physics. 

Meson and baryon spectra are well described in the constituent quark model (CQM). 
The CQM Hamiltonians written in the 1960's contained only the kinetic terms and short 
distance spin-spin interaction. They quite successfully predicted the magnetic moments of 
the ground state baryons, as long as the magnetic moments of constituent quarks had their 
classical values. The use of a nonrelativistic model was justified for cc and bb mesons, but 
did not work very well for light mesons. More sophisticated Hamiltonians treated spin- 
dependent interactions nonperturbatively, and based them on QCD [Hj. It was possible to 
describe hadrons within a unified, relativized quark model with chromodynamics, in which 
the qq interaction is a sum of the Coulomb (one-gluon-exchange) potential and a linear 
confining term expected from QCD 0IB|- 

High energy hadron-hadron and hadron-nucleus scattering at small and intermediate 
momentum transfers are well described by assuming that mesons and baryons are bound 
states of two and three constituent quarks, respectively [S]. Moreover, exclusive processes 
at small and intermediate momentum transfers agree well with the constituent quark model 
predictions of elastic and transition form factors. Therefore, the CQM approach provides 
a relevant description of nonperturbative QCD at low and intermediate momentum scales. 
Because it is hard to derive the constituent quark model from QCD, one may search for a 
relation between CQM and sum rules |10j . 

1.2 Mesons 

The most convenient formulation of QCD is a constituent representation in which hadron 
states are dominated by a small number of constituents. It will be assumed that in this 
representation, interactions that change the number of particles as well as other relativistic 
effects are small. A natural choice is the Coulomb gauge because it operates with simple 



degrees of freedom in the nonrelativistic limit. This framework works especially in QED, 
where for example the hydrogen atom is very well described by the Coulomb potential. The 
Coulomb gauge will be discussed more in Chapter 3. 

Each hadron is composed of quarks and may contain valence gluons. The simplest 
configuration of quarks that gives a color singlet is a quark-antiquark pair (meson). The 
next possibility for a colorless strongly interacting particle is a bound state of three quarks 
(baryon). One could construct more complicated configurations, for example the so-called 
pentaquarks, but until now, there has been no strong evidence for objects built of more 
than three quarks. 

Mesons can be classified based on their quantum numbers J PC . Here J is total angular 
momentum of a particle, P is its parity, and C denotes charge conjugation. The total 
angular momentum is given by 

J = L + S, (1.1) 

where L is relative orbital angular momentum of a quark-antiquark pair and S denotes total 
intrinsic spin of this pair, 

S = S! + S 2 . (1.2) 

Because quarks are fermions with spin 1/2, the values of S can be either or 1. Thus, the 
values of J are integer and mesons are bosons. The orbital angular momentum L can take 
any integer positive value or zero, and this determines all possible meson configurations. 

Parity, which determines how the sign of the wave function of a particle behaves under 
a spatial reflection, can be obtained from 

P = (-l) L +\ (1.3) 

whereas charge conjugation, describing the particle-antiparticle symmetry (well-defined only 
for neutral mesons composed of a quark and an antiquark of the same flavor), is given by 



C = (-l) 



(1.4) 
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Table 1.1: Classification of the simplest light imflavored mesons. 
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Table 1.2: Flavor SU(3) pseudoscalar and vector meson nonets. 

In the case of two-boson systems, parity would be given by a slightly different formula, 
P = (-l) L . 

The u flavor is given the third component of isospin /3 = 1/2, whereas the d flavor 
has I3 = —1/2. The concept of isospin came from the idea of treating the proton and the 
neutron as two states of one particle, the nucleon, having two values of the ^3, like fermions 
have two values of spin quantized along a fixed axis. Light unflavored mesons, i.e., mesons 
containing only flavors u and d, have I3 equal to either or 1. For example, the pion with 
J PC = ^ and 1=1 has three isospin components: 7r + , 7r° and 7r _ (isospin triplet), 
whereas for 1 = there is only one component: r] (singlet). The s quark has isospin 0; thus 
strange mesons have isospin 1/2. From three quark flavors we can build up nine mesons 
grouped into an octet and a singlet. In Table 11.11 we present the classification of light 
unflavored mesons with respect to the above quantum numbers. In Table lL2l we show the 
flavor wave functions of the nine pseudoscalar (J PC = '") and nine vector 

(J PC = 

mesons. 

In reality however, we do not observe exact configurations corresponding to t]q and rjg 



but rather their linear combinations known as rj and rf . The transformation matrix between 
both pairs must be orthogonal and thus has one parameter, the mixing angle 9. A similar 
situation occurs for loq and u)g, but in this case the mixing angle is such that the observed 
particles are given by 



For a particular combination of L and S we have more than one flavor multiplet due 
to different radial quantum numbers. Normal pions and kaons have n = 1, and are the 
lightest, whereas radially excited mesons with higher values of n are heavier. For example, 
the 7r(1300) is a candidate for a radially excited ir, and the /o(1450) is a candidate for a 
radially excited p. In this work we will deal only with radial ground state mesons because 
their orbital wave functions should have the same size as the well-known pions and kaons. 
Knowledge of a meson's structure is the first step towards understanding its dynamics, 
which is responsible for the spectrum of all observed mesons and their decay widths. In 
fact, all mesons that can decay strongly are not bound states but resonances, and we can 
study QCD by analyzing how they decay. The quark model description of such decays 
assumes qq pair creation in the gluonic field of the decaying meson. A phenomenological 
model based on quark-antiquark pair production from the vacuum is referred to as the 3 Po 
model ^21 El d| El El d| ■ Although this decay mechanism gives quite successful values 
for meson widths, it is not rigorously related to QCD which allows qq creation only from 
a gluon. The recent benchmark predictions for decay widths of light mesons are given in 




1 



(1.5) 
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Chapter 2 



Exotic mesons 



In the preceding chapter we presented classification of the low-lying unflavored mesons. 
There are, however, certain combinations of internal meson quantum numbers: spin J, 
parity P, and charge conjugation C, which are missing in this classification, such as 
, H , 1 or 2~~. These quantum numbers cannot be obtained from adding the 
quantum numbers of the quark and the antiquark alone. The corresponding mesons are 
referred to as the exotic mesons. 

From lattice QCD and model calculations it follows that the lightest exotic mesons may 
be obtained by adding an extra constituent gluon with J PC = 1~~ to a quark-antiquark 
system. Such qqg states are referred to as hybrid mesons. In this work we will focus on 
mesons with the J PC = 1 h quantum numbers. The isovector multiplet with J PC = 1 h 
71"+, 7pf, -Ki, is predicted to be the lightest exotic One must emphasize however, that 
qqg states can also have nonexotic quantum numbers. The hybrid components of normal 
mesons may be important in the mechanism of meson decays. Nonexotic hybrid mesons 
will be discussed in Chapter 7. 

Hadrons with excited gluonic degrees of freedom may supply new insight into quan- 
tum chromo dynamics at low energies, where the gluon dynamics should be responsible for 
phenomena such as color confinement and dynamical symmetry breaking. Therefore, the 
discovery of exotic mesons is of a great importance. In this chapter we will briefly review 
the experimental situation in the search for the 7T±, and present theoretical predictions for 
its mass and width. 



2.1 Experimental situation 



A resonance can be identified by analyzing the spectrum of its decay products. If a strong, 
narrow resonance is present, this dependence takes the form of a sharp peak, as shown in 
Fig. 12.11 In this picture presenting the BNL E852 data of the r]ir channel in the charge 
exchange reaction ir~p — » ry7r°n |29j . two resonances can be clearly seen. The corresponding 
particles are the ao(980) and the 02(1320). If a resonance is weakly produced, an amplitude 
analysis may be required which identifies the resonance by a phase motion of the amplitude 
as a function of the invariant mass of the decay products. 
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Figure 2.1: Distribution of the r/7r° effective mass in the reaction 7r p-> rjTt°n at 18 GeV. 
The two resonances correspond to the ao(980) and ei2(1320) mesons. 



Using such amplitude analyses, several candidates for the tti have been recently reported. 
The 7Ti(1400) with mass M = 1370 ± 16tj$ MeV and width T = 385 ± 40l^ 5 MeV, was 



reported by the E852 Collaboration in the ??7r~ channel of the process ir~p — > 7r~r]p |20|I21|. 
This state was confirmed by the Crystal Barrel Collaboration in the nr] channel in the 
reactions pn — > tt~tt°t] |22j and pp — > 7r°7r°77 |23j . In the r]Tr° channel, two resonances shown 
in Fig. 12.11 provide benchmarks for the amplitude analysis. A possible signal on the order 
of 1% of the dominant <22(1320) have been extracted in both rjit and r]Tr~ final states. 

A Breit-Wigner (BW) parametrization of the S-wave and D-wave corresponding to the 
ao and mesons in the 7/71"° and rjir~ channels is confirmed by the data, but the resonance 
interpretation of the P-wave is problematic. First, the left panel of Fig. l2~2l which represents 
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Figure 2.2: Comparison of the data on r]n~ (left) and t]'tt~ (right) production in the n~p 
interaction at 18 GeV with the results of the amplitude analysis: the D-wave (dashed line), 
the S-wave (dash-dotted line), and the exotic P-wave (dotted line). The resonance peaks in 
the D-wave correspond to the a2(1320). 

the t]tt~ spectrum shows that the signal for the 7Ti(1400) is weak. Second, it is impossible 
to find a selfconsistent set of the BW parameters for the P-wave. As a result, its phase 
as a function of the invariant mass does not increase over 90° which is required for a 
resonance |2*51 l2*§] . 

The E852 Collaboration has also reported two 7Ti(1600) states. One of these has M = 
1597 ± 10±fo MeV and r = 340 ± 40 ± 50 MeV, and decays into rfn 24 . In this channel, 
two strong amplitudes are extracted corresponding to the 02 (1320) and 7Ti(1600), as shown 
in the right panel of Fig. 12.21 The exotic signal is here much stronger, as compared to those 
in the rjn° and rjn~ channels. 

The other tti(1600) state with M = 1593 ± 8^ MeV and r = 168 ± 20±^|° MeV, 
was reported in the p°n~ channel |25| I26j . In this case, all expected well-known states: 
ai(1260), a2(1320), and 7^(1670) are observed, as shown in Fig. 12.31 25 . In addition, the 
amplitude analysis shows that the amplitude with exotic numbers J PC = 1 h has structure 
which is consistent with a resonance at 1.6 GeV decaying into pn. Evidence for the 7Ti(1600) 
has also been reported by the VES collaboration in three channels, bin, rj'n and pn |27j . 
with M = 1.61(2) GeV and T = 0.29(3) GeV. The tti(1600) signals in all these channels 
are somewhat different from one another, and therefore further experiments are needed to 




clarify the nature of these signals. 
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Figure 2.3: Distribution of the effective mass in the reaction ir p ^ ir + ir it p at 18 



Only the 7Ti(1600) reported in the pir channel has a width on the order of 100—200 
MeV, i.e., comparable to other meson resonances. The broad structures in the rjir and rfir 
channels can be accounted for by low-energy rescattering effects [3Uj. It is possible however, 
that the 1 ^ exotic meson in the rjiz channel is the same as 7Ti(1600) meson seen through 
its decay into pit. However, at this point this is only speculation [2*5] . 

2.2 Theoretical predictions 

The mass of the %\ can be obtained from calculations based on lattice QCD |391 1401 HT1 142j . 
They give values in the region 1.8 — 2.0 GeV. Theoretical predictions for this mass are based 
on various models. The QCD sum-rule predictions vary widely between 1.5 and 2.5 GeV 
ED EH EH]- The MIT bag model places this mass in the region 1.3-1.8 GeV ISB] . 
According to the constituent gluon model, light exotics should have masses in the 1.8 — 2.2 
GeV range [37j. The diquark cluster model predicts the tt\ state at 1.4 GeV (HE]- Finally, 
the flux tube model predicts the 1 ^ mass similar to the lattice results |491 1531 131] . 

In Tables 12.11 and 12.21 we show tt\ width predictions calculated using various nonrela- 
tivistic models: IKP CP [Mj and PSS JSHj. At 7ri mass equal to 1.6 GeV, the dominant 
modes are 7r&i(1235) and 7r/i(1285). For larger values of this mass, the above modes are 
still dominant, together with the i^^i(1400). 
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Table 2.1: Decay widths in MeV of various nonrelativistic models for the tti with mass 1.6 
GeV. 



m -> 


vr6i(1235) 
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Table 2.2: Decay widths in MeV of various nonrelativistic models for the tt\ with mass 2.0 
GeV. 

It should be noted that in each channel, one outgoing meson has orbital angular mo- 
mentum L = (these mesons such as n or K are called S-mesons). The other has L = 1 
(these mesons such as b\, f\ or K\ are called P-mesons). The above models favor modes 
that satisfy the so-called S + P selection rule. It states that a hybrid meson prefers to decay 
into one S-meson and one P-meson. There is a chance that relativistic corrections could 
significantly change this situation and favor the rjir, r/'n, and pir modes. This work aims to 
explore that possibility. 

Theoretical predictions indicate the importance of searching for the ~k\ in the p7r, bin, 
and /i7r channels. They also suggest a search for the KK\ channel. In order to compare 
these predictions with experiment however, more knowledge of branching ratios is necessary. 
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Chapter 3 



Dynamical foundations 



In the preceding chapter we described exotic mesons as quark-antiquark-gluon bound states. 
In order to proceed to their dynamics, we need to know how to obtain the exotic meson 
wave functions. In the nonrelativistic case, this can be done by using the Born-Oppenheimer 
approximation which works quite successfully for normal mesons treated as qq states. Fur- 
thermore, this procedure together with lattice simulations will provide some important 
information about the composition of the lightest hybrid meson. 

The constituent gluon plays a central role in the structure of an exotic meson. As the 
photon in QED, the gluon needs to be described in a particular gauge. A natural framework 
for introducing the constituent quark model and providing insights into calculating meson 
decays is the Coulomb gauge, which is free of unphysical degrees of freedom and has a good 
nonrelativistic quantum-mechanical limit. Relativization will be accomplished by using a 
relativistic phase space and transforming quark-antiquark states under Lorentz boosts. 

The model presented in this work is microscopic, i.e., at the level of quarks and gluons. 
However, mesons interact with each other via meson exchange. This force may contribute 
significantly to the dynamics of exotic mesons, and the quantitative analysis of this problem 
will be the subject of Chapter 9. 

We will begin the present chapter with the Born-Oppenheimer approximation. Then we 
will review the Coulomb gauge for QCD and describe the Coulomb gauge picture of normal 
and hybrid mesons. Finally we will proceed to relativistic effects in the Coulomb-gauge 
constituent quark model. 



3.1 Heavy quarkonia and the Born-Oppenheimer approxi- 
mation 



Lattice simulations are quite successful in predicting mass spectra for mesons and baryons. 
Hadronic decays, however, provide the real difficulty for such estimates. Thus we are left 
with phenomenological models of the coupling between mesons and hybrids. In general, 
there are two approaches for describing hadronic decays of hybrid mesons. The first regards 
a hybrid as a quark-antiquark state with an additional constituent gluon |4Sj . Such a 
meson would decay through gluon dissociation into a qq pair 03] . The second approach 
assumes that a hybrid is a quark-antiquark pair moving on an adiabatic surface generated 
by an excited gluonic flux-tube [310 02]. I* 1 this case a hybrid meson would decay because 
of phenomenological pair production described by the 3 Po model 03 OHl 021 021] • Recently, 
an extended version of the flux tube model has been introduced |58j . 

If constituent quarks composing mesons are heavy, such systems (heavy quarkonia) can 
be studied using the Born-Oppenheimer approximation [IB] E3 HE]- I n this approach it 
is assumed that formation of gluonic field distributions decouples from the dynamics of 
the slowly moving quarks, and therefore hadronic decays can be described within nonrela- 
tivistic quantum mechanics. This approximation can be justified for light quarks, because 
dynamical chiral symmetry breaking leads to massive consituent quarks. 

In the Born-Oppenheimer model, a hybrid meson is treated analogously to a diatomic 
molecule in which the heavy quarks correspond to the nuclei and the gluon field corresponds 
to the electrons. Initially, a quark and an antiquark are treated as spatially fixed color 
sources and this determines the glue energy levels as a function of the qq separation. Each 
energy level defines an adiabatic potential V q q(r). The quark motion is restored by solving 
the radial Schrodinger equation for each of these potentials. 

The lowest static potential gives a normal meson spectrum, whereas the excited poten- 
tials lead to hybrid mesons. The static potentials are determined from lattice simulations. 
The gluonic configurations can be classified according to symmetries of the qq "molecule". 
The strong interaction is invariant under rotations around the qq axis, a reflection in a plane 
containing the pair, and with respect to the product PC. Each configuration can be thus 



1 2 3 4 5 
r/r 

Figure 3.1: The static quark potential V s +(r) and selected gluonic excitations. 

labeled by the corresponding eigenvalues, denoted by A (the magnitude of the projection 
of the total gluon angular momentum onto the molecular axis), Y = ±1 (the sign of this 
projection), and PC = ±1, respectively. 

States with A = 0, 1, 2, ... are denoted by X, II, A, respectively. States which are even 
(odd) under the combined PC operation are denoted by g (u). Lattice simulations for the 
ground state configuration and the lowest gluonic excitations are shown in Fig. 13.11 |46| . 
The parameter ro is on the order of 0.5 fm. In the ground state (normal meson) A = 
(Eg ) and for the first excited state A = 1 (H u ). 

If the gluon is in a relative S-wave with respect to a qq pair, it has PC = +1. Lattice 
results show, however, that the lowest excited configuration has the gluon with PC = — 1 so 
the gluon orbital angular momentum with respect to a qq pair must be odd. The simplest 
choice is L = 1. Therefore, in Chapter 5, in order to construct the tt\ spin wave function 
we will couple a transverse gluon to the qq state with the p quantum numbers, in a relative 
P-wave. 

3.2 The Coulomb gauge 

In quantum electrodynamics, the electromagnetic field arises naturally from demanding an 
invariance of the action under the local gauge U(l) transformation. If (j> is a complex scalar 



field then the corresponding Lagrangian 



L = (S^)(^)-mV^ (3-1) 



is invariant under the transformation 



(ff^e-^fa cf>* -»• e*V, (3.2) 

where A is an arbitrary real constant. If A depends on the spacetime coordinates, however, 
then the derivative of the field does not transform covariantly, i.e., in the same way as (p. 
In order to remedy this problem one introduces the covariant derivative (like in the general 
theory of relativity) 

D„ = d lx + ieA lx , (3.3) 

where e is a real constant (the electric charge) and is the electromagnetic potential. 
This potential must transform according to 

A lt ^A ll + - e d^k. (3.4) 

The quantity 

= d fl A u - d v A fl (3.5) 

is the electromagnetic field tensor and its six nonzero components correspond to the fields 
E and B. The simplest Lagrangian built up from the gauge-invariant quantities is thus 



L - -\F^\ (3.6) 



and its variation with respect to A^ leads to the Maxwell equations in vacuum. Because 
of the gauge invariance we may introduce one constraint on the components of the field 
potential. In the Coulomb gauge this constraint is given by 

V • A = 0. (3.7) 
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In quantum chromo dynamics, the local gauge transformations form the SU(3) group 



<f>i Sirfj = (e5 A )^-, (3.8) 

where i, j = 1, 2, 3 and A; = 1..8. The matrices A are the hermitian and traceless generators 
of SU(3) (Gell-Mann matrices). In this case the expression for the covariant derivative is 
given by 

\ k 

D^dp + igA* — , (3.9) 
whereas the potential transforms according to 

4y - sa »y s ~' ~ l -(d,s)s-\ (3.io) 

The gauge-invariant field tensor is given by 

G%, = d,A a u - d„A% + gr hc A\Al, (3.11) 
where f a b c are the SU(3) structure constants, and the simplest field Lagrangian is thus 

L = -\g%G^. (3.12) 
The chromoelectric field corresponds to the Go a components of the field tensor 

and satisfies the Gauss law 

V -E a + gf abc A b -E c = gp a q . (3.14) 

Here = ifi^-ip is the quark color charge density. Introducing the covariant derivative in 
the adjoint representation 

Df = 6 a % + igT^ b Al, (3.15) 



where T Q c b = if cab , leads to 

D a6 • B b = gp a q . (3.16) 
If E|| = — V0 is the longitudinal part of of the chromoelectric field then we obtain 

- (D ab • V)/ = gp\ (3.17) 

where p = p q + p g is a total color charge density. Here the transverse gluon color charge 
density is given by 

where Ej_ is the transverse part of the chromoelectric field. Combining the above equations 
leads to: 



1 gp\ 



V D 



A 0a = tt^-V 2 )— i— gp a . (3.19) 
The last equation results in the instantaneous nonabelian Coulomb interaction Hamiltonian 

R c = \J d 3 xd 3 yp a (x)K ab (x,y;A)p b (y), (3.20) 



where 



K ab (x,y;A) = ( X , a \^(-V 2 )^\y,b). (3.21) 



After quantization, the field Ej_ becomes the momentum conjugate to the vector potential. 
The confining nonabelian Coulomb potential will be represented in the diagrams below by 
the dashed lines. More details related to properties of the Coulomb-gauge QCD can be 
found in Ref. [TT] . 

A simple phenomenological picture of hadrons and their decays in terms of quantum 
mechanical wave functions emerges naturally in a fixed gauge approach. In the Coulomb 
gauge, for example, the precursor of flux tube dynamics originates from the nonabelian 
Coulomb potential, which also determines the quark wave functions |1 11 157j . The string 



couples to a qq pair via transverse gluon emission and absorption and such a coupling carries 
the 3 S± quantum numbers. 

In a description of decays based on Coulomb gauge quantization it is necessary to include 
the hybrid quark- antiquark- gluon configurations, since they appear as intermediate states in 
the decay of mesons. If such hybrid states also exist as asymptotic states, they would provide 
insight into the dynamics of confined gluons [82] . Fig. 13.21 shows diagrams corresponding 
to strong decays of qqg hybrid mesons (top) and qq normal mesons (bottom). The gluons 
connecting the Coulomb lines represent formation of the flux tube, e.g. the gluon string in 
the ground state. The overall initial state is enclosed by the solid oval. In the lower diagram 
the hybrid meson state appears as an intermediate state in a normal meson decay, which 
is assumed to proceed via mixing of a qq pair with a virtual excitation of a gluonic string 
and its subsequent decay. 




Figure 3.2: Strong decay of a hybrid meson (top) and a normal meson (bottom). 

In the Coulomb gauge the quantum numbers of the gluonic states can be associated with 
those of a transverse gluon in the presence of the static qq source. This is because transverse 
gluons are dressed I74j. and on average behave like the constituent particles with the 
effective mass m g ~ 500 MeV [23 US]- Thus low-energy excited gluonic states are expected 
to have a small number of transverse gluons. The flux tube itself is expected to emerge from 
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Figure 3.3: A three-body potential between a transverse gluon and a static qq source. 

the strong coupling of transverse gluons to the Coulomb potential. The transverse gluon 
wave function can be obtained by diagonalizing the net quark-antiquark-gluon interactions 
shown in Fig. 13.31 (in addition to the gluon kinetic energy) . 

The above three-body interaction plays an essential role in the dynamics of hybrid 
mesons. A transverse gluon has a gradient coupling to the Coulomb potential. Thus the 
P-wave transverse gluon receives no energy shift from this coupling and the energy of the 
S-wave gluon state is increased. In the Coulomb gauge picture, the shift of the S-wave state 
via this three-body interaction may be the cause of the inversion of the S— P levels seen on 
the lattice. Using only two-body potentials between quarks and gluons leads to the S-wave 
gluon in the lowest energy excited state, which disagrees with lattice data |59| . 

3.3 Relativistic effects 

Another and very important issue is the question of relativistic effects. Even though a 
simple nonrelativistic description appears to be quite successful in predicting decay widths 
of mesons as heavy as ~ 1 — 2 GeV, the presence of light quarks raises the question of 
validity of this description. It has been shown that relativistic effects for hadronic form 
factors may be significant |6U1 1611 WR I63j . It is possible that they are responsible for the 
discrepancy between experimental data and theoretical predictions. This work will try to 
estimate the size of relativistic effects applied to the tt\ decays. 

In order to calculate relativistic decay amplitudes exactly in the Coulomb gauge, one 
needs to find the fundamental quantities (dynamical generators of the Poincare group) in 
terms of the chromodynamical fields. This problem, however, is very difficult to state and 



in this work we will not solve the Coulomb gauge QCD Hamiltonian to obtain meson wave 
functions. Instead, we will use the general transformation properties under the remaining 
kinematical symmetries (rotations and translations) to construct the states. 

The relativistic meson and hybrid spin wave functions will be elements of irreducible, 
unitary representations of the Poincare group for noninteracting particles. The interaction 
between particles should enter the dynamics by finding a new mass operator (the Bakamjian- 
Thomas model). Because an exact form of the strong potential between quarks is unknown, 
we will employ instead a simple parametrization of the meson orbital wave function. This 
is clearly an approximation which cannot be avoided without solving dynamical equations 
for the boost generators E2] . 



Chapter 4 



Relativistic dynamics 

In this chapter we will review the Lorentz group in vector and spinor representations, which 
is a framework for the ten fundamental quantities describing the dynamics of a system of 
noninteracting particles. Then we will introduce the Bakamjian-Thomas construction of 
these generators for interacting particles. Finally we will discuss a Wigner rotation, which 
plays an essential role in constructing relativistic and covariant spin wave functions for 
mesons and hybrids. 

4.1 The Lorentz group 

The principle of relativity in the framework of general relativity requires that physical 
laws must be invariant under all transformations of the coordinates. Gravitational fields 
are automatically included if one deals with curvilinear coordinates, however, they are 
important only for large-scale phenomena. Yet in the physics of elementary particles, the 
curvature of the spacetime is small and can be neglected. Therefore one needs to deal 
only with the metric tensor of a flat spacetime. In this case the principle of relativity 
requires that physical laws must be invariant under transformations from one inertial frame 
to another. Such transformations are called inhomogeneous Lorentz transformations, and 
the coordinates transform linearly according to 

x llx = A^x" + a", (4.1) 



where A^, is the Lorentz matrix, and is a constant four- vector. From the invariance of 



the finite interval X ^ X ^ — X^ X ^ for a^ 1 = 0, it follows that the Lorentz matrix must be 



orthogonal, 



(4.2) 



or AA T = 1. Thus its determinant can be either 1 or -1. 

All inhomogeneous Lorentz transformations can be divided into four categories, depend- 
ing on the signs of the determinant of A and the component A° . We will be interested in the 
proper Lorentz transformations, having both signs positive. They can be built up from in- 
finitesimal transformations involving boosts, rotations and translations, but cannot involve 
reflections. The proper inhomogeneous Lorentz transformations are continuous and form a 
Lie group called the Poincare group. If = then this group is called the Lorentz group. 
The principle of relativity will be satisfied if physical laws are invariant under infinitesimal 
transformations given by (|4.1j) in which 



(4.3) 



where lo^ v are infinitesimal quantities that are antisymmetric oj^ v = —uj u ^. This property 
results from the orthogonality of A^- 

Rotations are orthogonal transformations of the coordinates mixing their spatial com- 
ponents and form a subgroup 0(3) of the Lorentz group. They are described by a Lorentz 
matrix with A° = 1, A°j = and A* = 0. The remaining components are functions of 
three angles which may be chosen as the Eulerian angles of a rigid body. For example, a 
rotation by the angle 4> about the z-axis corresponds to the Lorentz matrix, 



R z (<t>) 



( i o 

cos 4> 

— sin i 

\0 





sin< 
cos i 







1 J 



(4.4) 



and similarly for two other axes. For an infinitesimal angle of rotation we can write only 



linear terms in (f), 



R z (5cp) = l + iJ z 5(f>, 



(4.5) 



and the passage to a finite rotation is given by 



R Z (<P) = lim {\ + iJ z ^-) N = e 



N _ JJ z <f> 



N^oo 



N J 



(4.6) 



The matrix J z is called the generator of the rotation about the z-axis. The rotation group is 
nonabelian, i.e. [R x ,R y ] ^ [R y ,R x ]. The corresponding generators satisfy the Lie algebra, 



[Jii Jk] — i^iklJl- 



(4.7) 



Boosts are described by a Lorentz matrix with A l j = 0. For example, a boost in the 
z-direction with velocity v corresponds to 



^ cosh-;/; sinh?/^ ^ 

10 

1 

t smh^/j cosh^ J 



(4- 



where cosh-0 = 7 = (1 — v 2 ) 1 > 2 . Writing 



(4. 



for an infinitesimal value of ip leads to the Lie algebra of the homogenous Lorentz group 
(rotations and boosts), given by (|4.7[) and 



[Ji,Kj] = ieij k K k , 
[Ki, Kj] = -ieij k J k . 



(4.10) 



o ,1 



An arbitrary vector (p°,p) transforms under the boost with velocity v according to 

p' = 7(p° + v-p), 



P|| = 7(P|| + V) 
P± = P±, (4.11) 

where pn + pj_ = p and pn = (v • p)v/u 2 . Parametrization by M, P and E such that 



V7 = P/M and E = y 'M 2 + P 2 leads to the following transformation laws: 



P'° = P"t7 + 



oE , P-P 



M M ' 

p' = p + p °P + <P P ' P V (4.12) 

If we introduce the four-dimensional antisymmetric generators defined by 

Jij = —UjkJk, 'hi = Ki, (4-13) 
then the Lie algebra may be written as one equation, 

\Jfiui Jpa] — i>\9vpJ ~\~ 9p,crJup 9p,pJucr dvoJpp)- (4-14) 

The components J^ u are proportional to the quantities co^ given in (|4.3jl . and infinitesimal 
constants of proportionality are either 5<j> or dip. 

The above operators may be expressed as differential operators instead of matrices. 
This will enable us to introduce the generators of translations. For an infinitesimal rotation 
about the z-axis we have 



T tl s r -fix ,y ,z)- f(x,y,z) 

Jzf{x,y,z) = bm? = 

<p— >o (p 

v ■ f(x + y<f>,y- x(f>, z) - f(x, y, z) d d 

= JST 4> = l{y Tx ~ X d~y )f - (415) 



r\ r 



For a boost in the z-direction we obtain similarly, 

■§;>■ <416) 

One can easily check that such defined operators satisfy the Lie algebra of the homogenous 
Lorentz group, (|4.7jl and I4,1U|) . For translations we can write 

rr fU \ r J(t,x,y,z + ()- f(t,x,y,z) .d 

T z f(t,x,y,z) = hmt = i—f, (4.17) 

and this leads to 



[T„T v ] = 0, 

[T p , Jfiu] = i{ T v9w ~ T ii9v P )- (4.18) 
Equations (|4.14l) and Q4.18JI constitute the complete Lie algebra of the Poincare group. 



4.2 The ten fundamental quantities 

Another requirement for a dynamical theory is that the equations of motion should be 
expressible in Hamiltonian form. This is necessary in order to make a transition from 
classical to quantum theory. The dynamics of a system is described by quantities called 
dynamical variables, which for particles can be taken as their coordinates and momenta, 
and for fields as their four-coordinates in spacetime. Any two dynamical variables £ and 
r\ must have a Poisson bracket [£,r/], and its form must not change under an infinitesimal 
Lorentz transformation. From this it follows that each dynamical variable £ will change 
according to 

£' = e+[e,n (4.i9) 

where F is an infinitesimal dynamical variable independent of £ and depends on the change 
in the coordinate system. Thus it must depend linearly on the infinitesimal quantities uj^ 



and a a . Therefore we can write 



F = -P% + -M^uv, (4.20) 
where P^ and M^ v = —M u ^ are finite dynamical variables called the fundamental quantities 

Each of the ten fundamental quantities is associated with an infinitesimal transfor- 
mation of the Poincare group. Pq is the total energy of the system and is related to a 
translation in time, Pi form the three-dimensional total momentum and are related to 
translations in space, and My correspond to the total angular momentum and are related 
to three-dimensional rotations. The quantities Moi correspond to boosts but do not form 
any additive constants of motion. From the commutation relations between infinitesimal 
transformations Q4.14JI and (|4.18|) . it follows that the Poincare algebra is given by: 

[P^Pv}=0, 

[P P , M^v] = g vp Pp - g^pPv, 

[M^, M pa ] = -g pp M va - g va M w + g vp M^ + g^M vp . (4.21) 

In order to describe a dynamical system one must find a solution of these equations, i.e. P p 
and M^y. This is the central issue in relativistic quantum mechanics. 
A simple solution of (|4.21j) for a single point particle is given by 

P fl =p fl , M^ u = q^p v = q u p^, (4.22) 

where q^ are the coordinates of a point in spacetime and p p are their conjugate momenta, 

[Qn,Qv] = 0, [Pp,Pu] = 0, [Pp,q u ] = gpw (4.23) 

One usually works with dynamical variables referring to a particular instant of time. The 
fundamental quantities associated with transformations that leave this instant invariant 
(spatial translations and rotations) appear to be simple, whereas the remaining Pq and Mio 



called Hamiltonians are not. Without loss of generality we may take qo = 0. Therefore po 
no longer has a meaning. But we can modify formulae Q4.22JI in order to eliminate po from 
them. Let us take 

P l i=V l i J r \^{p p p P - m 2 ), 

M^v = q^Pu - q u Pp + ^iLv{v p V P ~ m 2 ), (4.24) 
where m is a constant, with an appropriate choice of and \^ u . This leads to 

Pi=Pi, Mij = qiPj - qjpi, 

Po = \JpjPj + m 2 , M i0 = qi^JpjPj + m 2 . (4.25) 

These are the fundamental quantities for a particle with mass m in the so-called instant 
form of dynamics. There are two other forms: the point form and the front form, but the 
quantities appearing there are not as intuitive as in the instant form |65| . 

If for a single particle we replace qi by the operators i-^- and Mjo by the so-called 
velocity operators Vi = |(%-ff + H<li), where H = Pq, we will transit to quantum dynamics. 
In vector notation we can write 

P = p, M = qxp, 

H = yjm 2 +-p 2 , V = ^(qH + Hq). (4.26) 
For two noninteracting particles the ten operators are given by sums, 
P = pi + p 2 , M = qi x pi + q 2 x p 2 , 

H= V / mf+p2 + y/mj + pl V = qiymf + pf + q 2 ^m 2 + p 2 . (4.27) 

The expression 

M = VH 2 - P 2 (4.28) 
is the mass operator of the system viewed as a single entity, and commutes with all ten 



operators l|4.27j) . 

The last topic we will discuss in this section is related to the Casimir operators of the 
Lorentz group, i.e., the quantities that commute with all ten fundamental quantities 
and M^ u . Using the equations of the Poincare algebra (|4.21|) one can show that the only 
operators that have this property are 



The first Casimir operator is just equal to M 2 , and is related to the mass of a system viewed 
as a single entity, or to the mass of a particle. In the rest frame of a massive particle with 
mass m, the operator C2 behaves like the square of the angular momentum operator M 2 , 
and for spin s has 2s + 1 eigenvalues. For a massless particle, however, there exist only two 
eigenvalues ±s 

4.3 Spinor representation of the Lorentz group 

We defined the Lorentz group via the transformation properties of the coordinates x^. 
Quantities that transform under the Lorentz tranformations in the same way as the coor- 
dinates are called vectors, and the matrix A^, is referred to as the vector representation 
of the Lorentz group. This representation is suitable when dealing with vector particles 
having integer values of spin. However, for particles with spin 1/2 (fermions) it is much 
more useful to introduce the spinor representation of the Lorentz group |73j . This will be 
the subject of the present section. 

Consider the group SU(2), consisting of 2 x 2 unitary matrices with unit determinant. 



Ci = P"i> c 2 = W^W, 



(4.29) 



where is the Pauli-Lubanski pseudovector 




(4.30) 



an 



These conditions imply 



U 



-b* a* 



, |a| 2 + |6| 2 = l. 



(4.31) 



It can be shown that if a matrix H is hermitian and traceless, so is H' obtained by the 
transformation H' = UHU^ . Let us introduce the matrix X given by 



\ 



x — ly 
y x + iy -z J 



(4.32) 



where cr,- are the Pauli matrices. Since X is hermitian and traceless, so is X' = UXUl We 
also have detX' = detX which gives 



12 , 12 i 12 2,2,2 

x +y + z = x + y + z . 



(4.33) 



This is the condition for a rotation of the position vector r. Therefore, we arrive at the 
conclusion that the SU(2) transformation is related to the 0(3) rotation. 

We would like to find the explicit form of the matrix U that corresponds to an arbitrary 
rotation. For a rotation about the z-axis we have 



x = x cos 4> + y sin (J), y = — x sin <p + y cos (J), z = z, (4.34) 



and substituting this into X'U = UX gives 6 = and a = e 1 ^! 2 . Thus 



e i<P/2 o 
e-^/ 2 



(4.35) 



This result can be generalized to a rotation about the axis with the unit vector n, 



U n {<P) = e^^ 2 = cos 4>/2 + i sin 0/2 a ■ n. 



(4.36) 



The above relation is very similar to the corresponding expression for the Lorentz matrix 



on 



for a rotation, R n ((j)) = e lJn , and this is related to the fact that the Pauli matrices satisfy 
the same commutation relations as the matrices Jf. 

[a i ,a j ] = ie ijk a k . (4.37) 

When a vector rotates by the full angle 2ir, a spinor rotates only by the angle it and changes 
sign with respect to the original value. Thus both matrices U and —U correspond to the 
same rotation matrix R. 

The matrix U is regarded as the transformation matrix of a two-dimensional complex 



object £ 



2 . 



£' = ft = j/tft. (4.38) 



The quantities having the above transformation property are called spinors. We see that 
f and f * transform in different ways, but we may show that and i^S, = 

transform in the same way under SU(2). We also notice that £^<j 2 f is a scalar under 
rotations, whereas £(i<J2&^ transforms like a vector. 

Now we proceed to transformations of spinors under boosts. From the Lie algebra of 
the Lorentz group (|4,7|) and Q4.1UJI it follows that the matrices K{ = ±iJi are its solutions. 
Therefore, spinors should transform under boosts according to formula (|4.36l) with the 
replacement o % — > ±ia l . We may define two types of spinors £ and rj, transforming with a 
plus or a minus sign, respectively. For the first one we have 

3^= a/2, K.W*>=-i*/2, (4.39) 

and if (<p, ip) are the parameters of a pure rotation and a pure boost this spinor transforms 
according to 

g = eur/Ht-i^t = C£. (4.40) 

For the second one we have 

jW=a/2, K^=ia/2, (4.41) 



o 1 



and the Lorentz transformation is given by 



These are inequivalent representations of the Lorentz group and there is no matrix S such 
that D = SCS" 1 . Instead, we have D = CT2CV2. The matrices C and D are no longer 
unitary, but still unimodular. Such matrices build the group SL(2,C) which is related to 
the Lorentz group like SU(2) was related to the rotation group. The matrix X is now given 
by 

X = 3%, (4.43) 
where o"o is the 2x2 unit matrix, and the transformation law has the form 

X' = GXG~\ (4.44) 

where G belongs to the SL(2,C). 

If we define the parity transformation v — ► — v, which changes the sign of K but leaves 
the sign of J, then the spinors £ and r\ will interchange. Therefore we may define the 
four-spinor ^ = , transforming under rotations and boosts according to 



e i/2a-(4>-iip) Q 

= | | (4.45) 

e i/2o--0+*V>) 



and under parity like 



U 1 



1 



(4.46) 



The bispinor ^ is an irreducible representation of the Lorentz group extended by parity. 
For pure boosts we can write 

£' = [coshV'/2 + a ■ nsinh^/2]£, 

rj = [cosh -0/2 - cr • nsinh^/2]?7, (4.47) 



where n is a unit vector in the direction of the boost. If £ and 77 refer to a particle at rest, 
then £' = £(p) and 7/ = r/(p), where cosh?/; = 7 and p = mv). One can show that in a 
moving frame of reference these spinors satisfy the equation 

/ 



— m E + a • p 
y E — a ■ p —777, 



*(p) = 0, (4.48) 
where E = \/m 2 + p 2 . Introducing the 4x4 Dirac matrices in the chiral representation, 

I • ~ ! ■ (4-49) 



(0 A 




( 















leads to the Dirac equation 

(7% - m)*(p) = 0. 
The matrices 7^ satisfy the anticommutation relation 

'f'f + VV = 2^, 



(4.50) 



(4.51) 



which is actually their definition. 

We will work in the standard representation of the Dirac matrices, in which 7 is diag- 
onal, 

/ 

1 
-1 



7 



(4.52) 



V 



It can be obtained from the chiral representation by 

7° = T 7 ° c T-\ 



(4.53) 



where 



T = 



V2 



1 1 
1 -1 



(4.54) 



00 



Therefore the bispinor becomes 



(4.55) 



and the spinor representation of the boost is given by the matrix 



cosh tb/2 sinh ib/2(j ■ n 

s-\ |. (i. on) 

smhib/2a • n cosh^/2 



or finally 



i / E(m, p)+m a ■ p \ 

5(0 -p)=5(m,p) = -====_ . (4.57) 

y/2m{E(m 1 p) + m) \ cr . p E(m,p)+mJ 

The spinor representation for rotations and boosts can be also derived from the assump- 
tion that the Dirac equation in the position space 

(i^dp - m)# = (4.58) 

is invariant under Lorentz transformations x^ — > x'^. For the bispinor we will assume that 
it transforms according to 

= 5* , (4.59) 

where S is a unimodular matrix corresponding to the Lorentz transformation. Its form can 
be derived from the requirement that the Dirac equation will remain unchanged, 

(i^'d^ -m)^' = 0, (4.60) 

leading to 

ff-VS = A^j u . (4.61) 

In order to derive the form of S, we will first consider the infinitesimal Lorentz trans- 
formation ()4.3|) . Consequently, the matrix S will be a linear function of the generators Uu V , 



and the simplest guess is 

S = l + auvYV, (4.62) 

where a is some constant. The inverse matrix S" -1 in the linear approximation which suffices 
for our considerations is 

S- 1 = 1 - ac^T'V- (4-63) 

Substitution of S and S^ 1 into (|4.58|) gives a = 1/4. The passage to finite transformations 
is again done by exponentiation, 

One can show that this expression is equivalent to (|4.57l) . 

The explicit form of bispinors in the frame of reference in which a particle has momentum 
p can be obtained by acting with (|4.57|) on the solutions of the Dirac equations in the rest 
frame, given by 

u(X)=( X(X) ), v(X)=( ° V (4.65) 

where A = ±1/2 are the values of spin. The term ia 2 in v guarantees that spinors v and v) 
transform in the same way under SU(2), and 

X T (+1/2) = (1,0), x T ("l/2) = (0,1). (4.66) 

Thus we get 

1 f(E(m,p)+m)x(X)\ 

u(p,A) = ^=== (4.67) 
y/E(m,p) + m V (a-p)x(A) / 

and 

1 / (o-"P)^2X(A) \ 

v(p, A) = — . (4.68) 
y/E{m, p)+m\ (E(m, p) + m)ia 2 xW / 

The nonrelativistic limit of the Dirac equation is free of paradoxical properties only in 
the first approximation. It is possible, however, to find a representation in which it is clear 
how to associate operators with classical dynamical variables so that these operators tend to 
their expected nonrelativistic form |7Uj . In the presence of an external field the nonrelativis- 
tic reduction is most conveniently obtained by an infinite set of canonical transformations 



related to a free field transformation. This is referred to as the Foldy-Wouthuysen repre- 
sentation, and can also be extended to Klein-Gordon and Proca particles In this work 
the nonrelativistic limit will be reduced to the first approximation. Therefore, the Dirac 
representation will be suitable for our purposes. 

4.4 The Bakamjian-Thomas model for interacting particles 

For a system with a fixed number of particles, Pj and Mij will be sums of their values for 
separate particles, 

Pi = J2 Pi, Mi^Y^iliPj-QjPi)- (4-69) 
For the Hamiltonians one must add the interaction terms, 

p o = Yl \JP3P3 + m2 + U > 

M i0 = U\fPjPj + m 2 + Ui. (4.70) 

From the commutation relations Q4.21JI it follows that U is a three-dimensional scalar, U is 
a three-dimensional vector, and 

Ui = qiU + bi, (4.71) 

where 6, is a constant three-dimensional vector. The remaining conditions for U and U 
are quadratic and therefore a construction of a complete dynamical theory of a relativistic 
theory is very difficult. 

An interaction enters only in H and V. A practical method of constructing the gener- 
ators in this case was developed in (HB], where the set of new operators satisfying simpler 
commutation relations was introduced. In this set, the interaction appears only in the mass 
operator (|4.28|) . Suppose we can make a transformation from qi, q2 and pi, P2 to the total 
momentum P, the coordinates of the center-of-mass R, the relative momentum p and the 
relative coordinate vector r. The commutation relations are not disturbed if: 

1. M = RxP + rxp, 

2. M depends on p only, 



3. V can be expressed in terms of M, P, M and R only. 

The interaction will be included if we replace M by any other function of p and r which 
is a scalar for space rotations. The only nonzero commutators of the set M, P, M, R are 



and the mass operator M is Poincare invariant if it commutes with P, M, R. Therefore 
it is only necessary to make sure that the above condition is satisfied. The macroscopic 
Hamiltonian of a system is given by 



and is obtained from the microscopic one (|4.7()|) via introducing the above relative variables. 
The above results can be generalized to systems with more than two particles, and to parti- 
cles with intrinsic spin |66| I68j . An explicit construction for a unitary operator that insures 
the free motion of the center of mass of any system is given in [HZj- However, for a given 
potential there is no unique way in which the relative variables may be defined _69; • Unfor- 
tunately, an exact form of the potential is not known and one must use phenomenological 
forms of the mass operator. In this work we will assume a gaussian form of the orbital 
wave function and fit it to a few measured form factors. This approach will not allow for 
a deeper understanding of the relativistic quark dynamics, although it makes possible to 
estimate relativistic corrections to the tt\ decay widths. 

4.5 Wigner rotation 

In this section we will derive how spin transforms under the boost transformations. It will 
be necessary for a construction of covariant spin wave functions for quark-antiquark pairs 
(mesons). Our goal is to solve 




(4.72) 




(4.73) 




(4.74) 



where A(p — ► q) boosts a particle with a momentum p to a frame in which the momentum 
is equal to q, and |p, A) is the spinor state. Multiplying the above expression by the identity 
A(0 — > p')A(p' — ► 0), where p' is obtained from p according to (|4.12|) . leads to 

A(0 - P)|p, A) = A(0 -> p')R(p, P)|0, A). (4.75) 

The quantity R is given by 

R(p, P) = A(p' - 0)A(p - p')A(0 -> P ), (4.76) 

and we will find its explicit form. 

In spinor representation we can write 

R(p, P) = S(m, -p')S(M, P)S(m, p), (4.77) 

where M parametrizes the boost A(0 — > P) like before. Substituting the expression Q4.57JI 
in the above equation leads after somewhat lengthy calculations to 

/ DW2)( Pi P) o 
i*(p,P)= |, (4.78) 

\ DWfaP) 

where 



< !, (q.P) 



(E(m, q) + m){E(M, P) + M) + P • q + ia ■ (P x q) 



y/2(E(m, q) + m)(E(M, P) + M)(E(m, q)£(M, P) + P • q + ml) 



AA' 

(4.79) 



It can be shown that L^ 1 / 2 ) has the form cos 0/2 + zsin</>/2<7 • n and thus the matrix 
H4.78JI represents a pure rotation. The matrix Z^ 1 / 2 ) is called the Wigner rotation matrix. 
In nonrelativistic quantum mechanics, spin should not change under boosts and this is 
reflected in the large-mass limit of formula (|4.79|) . 

D^^S XX , (4.80) 



O i") 



Finally, we obtain the transformation law for spinor states under boosts, 



A(0 - P)| P , A) = A(0 - p')D^\p,P)\Q,X') = D^(p,P)|A(0 -» P)p,A'>. (4.81) 

The state of a system having more than one spin index transforms like a spin tensor, i.e., 
each index transforms independently with the ZK 1 / 2 ) matrix according to formula (|4.81|) . 
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Chapter 5 

Relativistic spin wave function for 
mesons and hybrids 

Having described transformation laws for a single particle with spin 1/2, we may proceed 
to systems of noninteracting particles. We will focus on quark-antiquark pairs, i.e. mesons. 
This is necessary in order to construct the spin wave functions for the outgoing mesons 
resulting from the decay of the tv\. Since these mesons have nonzero momenta, a relativistic 
model of hadronic decays will have to include a Wigner rotation of spin. A similar con- 
struction is also required for the m spin wave function because the qq pair must be boosted 
to a moving frame before it couples with a gluon to a rest-frame hybrid. In the following 
section we will show how to build a relativistic spin wave function for each light unflavored 
meson (or a meson with equal masses of quarks) . Following that we will add the gluon and 
build the tt\. Mesons with different masses of quarks will be considered later. In each case 
we will begin with a rest-frame function, and then use the results of the preceding chapter 
to obtain a general expression for any frame of reference^. 

5.1 Meson spin wave functions 

The spin wave function for a meson is constructed as an element of an irreducible represen- 
tation of the Poincare group |fil | l72|. In the rest frame of a meson, the quark momenta are 
^This chapter is based on work by A.P.Szczepaniak 



given by 



Z£ = (£(m g ,q),q), = (E(m g , -q), -q), 



(5.1) 



and the normalized spin-0 and spin-1 wave function corresponding to J = + and 
are simply given by Clebsch-Gordan coefficients, 



(5.2) 



and 



*Jf = 0,o- g ,a 5 ) = (-,a g ;-,ag|l,A 9 g) = xV?)^|^xKK(A M -). (5.3) 

A factor io"2 accounts that the antiparticle spin doublet transforms under SU(2) in the same 
way as the particle doublet. The canonical polarization vectors 
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, e(0) = 





V o J 







(5.4) 



correspond to spin 1 quantized along the z-axis and satisfy the orthogonality relation 

£Y(Ay'*(A) = 5 ij . (5.5) 

A 

The invariant mass of the qq pair is 

rriqq = E(m q ,q) + E(rriq,—q), (5.6) 



where E(m, p) = ^/w 2 + P 2 ) and the total momentum of this system l qq is of course equal 
to zero. In the following we will assume 



Tn n = nia = m. 



(5.7) 



This condition is satisfied to a good approximation by the quarks u, d and may be used for 
a construction of light unflavored meson states. 

The rest frame wave functions (|5,2j) and Q5.3|) may also be expressed in terms of Dirac 



spinors quantized along the z-axis, 



V2m q g 



-u(q, <7 g )7 5 v(-q, cr 5 ), 



(5.8) 



and 



m. 



-u(q, cr g ) 



99 



2g* 



m 9 g + 2m 



uC-q^KCL^O.A-A (5.9) 



A spin wave function written in this form is manifestly covariant and thus it is straightfor- 
ward to find how it transforms under Lorentz transformations. In the above e l (l qq = 0, \ qq ) 
are the spatial components of the polarization four-vector 



e»(l qq = 0,\ qq ) = (0,e(\ qq )), 



(5.10) 



whose time component is zero in order to satisfy the transversity condition e^(k, X)ku = 0. 
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Figure 5.1: Boost of a quark-antiquark pair from the rest frame to a moving frame. 

Now we apply a boost from the rest frame of a qq pair to a frame of reference in which 
the momenta of the quark and the antiquark are \ q and l q , respectively and the total 
momentum is \ qq = L + l q , as shown in Fig. 15.11 The new momenta are given by 



ig = q + 



(q • W) 1 ^ 



E(m qq , Iqg) [m qq + E(m qq , \ qq ) 



m 



qq, 



(5.11) 



and 

!_ = _ q (SLiiilk i gM t (5.12) 

The spin wave function of a meson in a moving frame is obtained from the rest frame wave 
function, as we stated at the end of Chapter 4, by acting with the Wigner rotation matrix 
(|4.79j) on each of both spin indices. This gives 

9% (q,l 93 -,A ? ,A 9 -) = £ *Jf (q,V, = 0,a 9 ,a g -)D*^(q,l gg -)^(-q,l ?9 -). (5.13) 

The above Wigner rotation matrix corresponds to the boost with V7 = P/M. In spinor 
representation Eq. (|4.81j) leads to the following transformation laws for spinors vJ and v: 



T 9 



E ^ifCq. W)« f (q = ^ A 9 )5t(i 39 - - 0), (5.i4) 

where <S'(l g g — > 0) is the Dirac representation of the boost taking l q to q and l q to — q, given 
by ([4.57ft with p = — \ qq and M = m qq . From these laws we obtain the general form of the 
spin-0 wave function: 

9 qq (q, l qq , X q , X q ) = -J— U \l q , X q )S\l qq -> 0) 7 °7 5 5(1 OT - - 0)v(l q , X q ) = 

VZrUqq 

1 



u(l q ,X q )S-\l qq -> 0) 7 5 5(l 9q - - 0) W (1 ? , A,-) 



= - J— A g )7 5 f A 9 -). (5.15) 
Similarly we derive the spin-1 wave function: 



9 qq (q> Ag, Ag) 

V zm 9 g v rriqq + ±rn ' 

= -^—u(l q ,X q )S-\l qq - 0)( 7 4 - ^ )5(^ - 0)«(1 ? -, A 9 K(A M -) = 



1 



«(!,, A^Ayi™ - 0)7^ ^— Wl*, A 5 K(A OT -) 



99 ""ot 



m„o + 2m 



1 > 1 



-u(l q ,\ qJ 



qq 
1 



y tit -til 



y/2; 



m, 



-U(lq, X q ) 



VI 



m q q + 2m 



v(l q ,X q )A i 1/ (l qq ^O)e\X 



qq) 



q "q 



m qq + 2m 



v 0-q, X q )e^(l q q, X qq ), 



(5.16) 



where e M (l g g, X q q) are obtained from ()5.4I) through the boost with /?7 = \ qq /m qq : 



e 0-qqjXqq) 



Iqq ' e {Xqqj 



in 



qq 



e (l<?9> Xqg) — e(X q q) + 



(Iqq ■ e(Xgq))lqq 
m q q(E(m q q, Igq) + m qq ) 



(5.17) 



The invariant mass of the qq system is now 



mqq = m qq (\q, \g) = J (E(m, lg) + E(m, lg)) 2 - (l g + lq) 2 . 



(5.18) 



The wave functions (|5.15j) and IJ5.1fi|) are still normalized: 



Yj Vqqfa V x q)^qq(^ l qq> V x q) = *> 



£ *£"(q, W> A g , A $ )* g f (q, A„ A g ) = <^ A , ? . 



(5.19) 



By coupling the spin wave function (|5.8[) or (|5.9[) . respectively, with one unit of the 
orbital angular momentum L = 1, one obtains the rest frame spin wave functions for the 
quark-antiquark pair with quantum numbers J PC = l" 1 or ++ , 1 ++ and 2 ++ . Explicitly 
we have 



<f(q,W = o,^,a 9 -) 



N/2m g? (q, -q) 



«(q, a (? )7 5 ^(-q, cx 9 -)y 1Aq5 (q) (5.20) 



for the 1 + meson, and 



*Jf(q,W = 0i^»o-?) : 

xY u (q)(l,X;l,l\J,Xqq) 



2q i 



m q q + 2m 



o(-q,ff,)e < (0 I A) 



(5.21) 



for the J ++ (J = 0, 1,2). Here 1la( ( i) is a spherical harmonic and q = q/|q|. Using ()5.13j) 



i 1 



one can show that the wave functions for the qq pair moving with the total momentum l q g 
are given by 

tfjf l q , \ q , X q ) = 1 u(l q , X q h 5 v(k, A g -)riA M -(q), (5.22) 
V2m g q(l g ,l 5 ) 

and 

\J2m qq 1 m q q + lmi 

xy lz (q)(l,A;l,Z|J,A OT ->, (5.23) 

respectively, where m ? g = m 9 q(l 9 , lq) and q remains the same but must be written in terms 
of new variables: 



E( m, \ q )lqq (lqg • l q )l qq 

mqq m qq (E(m qq ,-\ qq ) + m qq ) ' 



q = A (l g( j - 0)1, = 1 9 - + ; - (5.24) 



In order to construct meson spin wave functions for higher orbital angular momenta 
L one need only to replace Y\\ with in ()5.22() and (|5.23() . In L = [formulae (|5.8j) . 
()5.9|) . ()5.15() and (|5.16[) ] we skipped a constant factor loo to make the spin wave function 
normalized to 1. But from now on, for consistency, we will assume this constant being 
implicitly included. In the nonrelativistic limit, where Wigner rotations may be ignored, all 
the spin wave functions for mesons simply reduce to the Clebsch-Gordan coefficients that 
we started from, coupled to appropriate spherical harmonics. 



5.2 The 7Ti spin wave function 

As we showed in Chapter 3, in the lightest hybrid meson wi wave function, a constituent 
gluon is expected to have one unit of orbital angular momentum with respect to a qq pair. 
Thus, the quantum numbers P, C require a quark and an antiquark to have parallel spins 
(5 = 1) 

In the rest frame of a 3-body system with a qq pair moving with momentum — Q and a 
transverse gluon with momentum Q, the total spin wave function of the hybrid is obtained 



by coupling the qq spin-1 wave function (|5.16|) and the gluon wave function ( J PC = 1 ) to 
a total spin S = 0, 1, 2 and J PC = 0++, 1++, 2++ states, respectively. Here, we will derive 
the expressions for each value of S separately, although the physical tt\ state should be a 
superposition of all three components. The way of calculating the corresponding coefficients 
in this linear combination will be given in Chapter 7. The total J PC = 1 ^ exotic meson 
wave function is then obtained by adding one unit of orbital angular momentum between 
the gluon and the qq: 

Aqg,Cr=±l,M,« 

xy„(Q)(5,M;l ) /|l,A ea: ). (5.25) 
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Figure 5.2: A constituent picture of a hybrid: the quark-antiquark pair quantized along 
the z-axis and the transverse gluon quantized in the helicity basis. 

The spin-1 rotation matrix relates the transverse gluon states in the helicity basis a 
(i.e., along its momentum) to the basis described by spin X g quantized along a fixed z-axis, 
as shown in Fig. 15.21 Thus, the helicity basis is rotated so that a coupling between a gluon 
and a qq pair may be done in the same basis with Clebsch-Gordan coefficients, since we are 
in the rest frame of a hybrid. Explicitly we have 



|Q,A 9 ) = ^Di 1) ;(0,0,-^)|Q,cT 



(5.26) 



where 9 and <f> are the polar angle and the azimuth of the direction of the gluon momentum 
Q, as shown in Fig. 15.31 For the gluon polarization vector with spin quantized along the 




Figure 5.3: A polar angle and an azimuth of the gluon momentum vector. 



z-axis we can write 



<7=±1 



where the helicity polarization vectors are given by 



(5.27) 



4(Q,a) = Y.D [ t ) M,-ct>yi\ g ) 



(5.28) 



Using the unitarity of the matrix D^> one can show 



(5.29) 



and with the help of the identity e^(Q, cr)e^(Q, a) = 5 13 — Q l Q 3 we finally obtain 



e*(Q,A 9 )=e^(A s )(^- QV 



(5.30) 



where Q i = QVIQI- 

The Clebsch-Gordan coefficients and the spherical harmonic in ()5.25j) can be expressed 
in terms of the polarization vectors (|5,4j) . For example: 



(1,A';0,0|1,A) =e*(A')-e(A), 
(l,A';l,A|0,0) = -i=e*(A / )-e*(A), 

(l,A';l,A"|l,A) = -^[e*(AOxe*(A")]- e (A), 



^(Q) = V^(0-Q- 



(5.31) 



,1 T 



Therefore we obtain: 

£<1, Xqq', 1, A 9 |0, 0)y u (Q)(0, 0; 1, l\l, X ex ) oc [e*(A^-) • e*(A g )][Q • e(A ex )], 
^(l,A^;l,A fl |l, S )y u (Q)(l,5;l,Z|l,A e:c ) oc [e* (Xqq) x e*(A g )] • [Q x e(X ex )}, 

l,s 

£<1, A OT -; 1, A fl |2, s)Yi,(Q)(2, s; 1, Z|l, A e:c ) oc Q • [e*(A 99 -) ® e*(A g )] • e(A« e ), (5.32) 

and the action of the rotation matrix on the gluon states results in replacing e l (X g ) 
with e* (Q, X g ). The normalized hybrid wave functions are then given by: 



qqg(S=0) y 8lT 

3 



£ <f (q, l qq = -Q, A ? , A,-)[e*(A w -) • e*(Q, X g )] [Q • e(A e x)], 



*fe(5=i) = V^E^f^ 1 ^ = -Q,A 9 , A 5 )[6*(A W -) x e*(Q, A fl )] • [Q x e(X ex )], 



qqg(S=2) y 
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^^f(q,l ?? - = -Q,A„A g -) 



xQ-[e*(A ?? -)®e*(Q,A ff )]-e(A e:c ), 



(5.33) 



where 



(A ® = i4iBj + AjBt - -6ij(A ■ B). 



(5.34) 



Writing the qq spin wave function more explicitly in terms of the quark momenta p,j 
and Pq gives 



^qTg(S) ( p <?' p «' "V X g) - - ^ 

X ll>n(S)(-Pq - Pq, A fl , A ex ), 



1 



i(Pg> \ 



P^-Pq 



VTlqq + 2m 



; (P?A) 



(5.35) 



where the gluon terms are respectively: 



^(s=o)(Q) A 9 , A ez ) = 
VV(s=i)(Q) A s , A e:r ) = 

■0/^(5=2) (Q, A 9 , A e:r ) = 



A e ^(Q,A 9 )Q'6'(A e:c ), 



M m q q(Eqq + TXl 



qq) 



e* c \Q,X g )Q«e l (X 



kji 



=i)(Q, A 9 , A e:r ) - -Vys=o) (Q> A fl , Aex)) , (5.36) 



and 



m q q = m q g(p q ,p q ), E qq - = E(m gq , Q), 

K = -Q, K° = E gg - + m qq -, e°(Q,A 9 ) = 0. (5.37) 

The loss of linear independence between all three functions Yws) came from the replacement 
of e l (Xg) with e*(Q, X g ), which is perpendicular to the momentum vector Q. 

The spin wave functions for other hybrid mesons can be constructed in similar fashion. 
In particular, in Chapter 7 we will describe the p and b\ mesons as gluonic bound states, 
and construct corresponding wave functions. 



A r\ 



Chapter 6 



Meson and hybrid states 

In the preceding chapter we constructed the spin wave functions for normal and hybrid 
mesons assuming that quarks do not interact, so each spin index can transform separately 
under a Wigner rotation. The interaction between a quark and an antiquark enters through 
the Hamiltonian H = P° and the boost generators of the Poincare group M . It is 
possible to produce models of interaction for a fixed number of constituents that preserve 
the Poincare algebra for noninteracting particles following the prescription of Bakamjian 
and Thomas, as we discussed in Chapter 4. Unfortunately, such a construction does not 
guarantee that physical observables such as current matrix elements or decay amplitudes 
will be relativistically covariant. Thus we must deal with phenomenological models of the 
quark dynamics, and in this chapter we will follow the common practice of employing a 
simple parametrization of the orbital wave functions. 

6.1 Mesons as qq bound states 

Unitary representations of noncompact groups are infinite-dimensional [M]. The rotation 
group is compact, because rotating by the angle 2tt (or 4ir for spinors) returns the trans- 
formed quantity back to the original state. However, this is not the case for boosts and 
therefore they do not form a compact group. This is reflected in the fact that the spinor 
representation of the Lorentz group (|4.45|) is not unitary. 

In quantum mechanics we are only interested in a unitary representation of a symmetry 



group, because the transition probabilities between states do not depend on the choice 
of a frame of reference. The problem of non-unitarity of the Lorentz group is solved by 
introducing the Fock space in which states are described by kets |p, A) with momentum 
p and spin A. This representation is infinite-dimensional because the spectrum of values 
of p is continuous, and thus it is unitary. It is also irreducible, because the states have 
well-defined values of mass m and spin s. In this section we will construct states for all 
mesons whose spin wave functions we have built in Chapter 5. 

The n (I = 1) and 77 (/ = 0) states {J PC = characterized by momentum P and 
spin Xgg, are constructed in terms of the annihilation and creation operators: 

io-+(p,/,/ 3)) = £ / {2 J*L p. ^nEL P ,) mm - p ' ) + E(m ' p?)) 



X (2^) 3 5 3 ( Pg + P g ~ P) — S c<iC9 (-,f q ; -,f q \l,I 3 )^ q g( Pq ,P q ,X q ,Xg) 

where the operators satisfy the anticommutation relations: 

{ & pA/c>kp'A'/'c'} = {V/c'VV/v} = (2vr) 3 2£;(m,p)5 3 (p - p )5 X \>5ff5 cc >, 
{bpXfcbp'X'fc'} = {dpxjcdp'x'fc'} = { b p\fci b p>\>f>ci} = { d P A/c d p>\>f>c'} = °> 
{bp\fc,dp>xf> c >} = {bp\fc,d p , x , flc ,} = {b p xf c ,dp'\>f>c>} = {bp X f c : d^yfc'} = 0.(6.2) 

In the above, ^f q g represents the spin-0 wave function (|5.15jl . written explicitly in terms of 
the momenta p q and Pg instead of the relativistic relative momentum q and the center-of- 
mass momentum P = \ qq (Eq. (|5. with \ q = p q and Eq. (|5. 12|) with 1^ = p^). The third 
component of isospin, flavor and color are respectively denoted by I3, f and c. The factor 
5 CqCjj guarantees that the meson state is colorless. 

The orbital wave function ipL results from the strong and electroweak interaction be- 
tween quarks that leads to a bound state (meson). Such a function depends on momenta 
only through the invariant mass of a quark-antiquark pair ()5.18|) . Normalization constants 
are denoted by ./V (with P = |P|) and the /i's are free parameters, being scalar functions of 



meson quantum numbers. Finally, the isospin Clebsch-Gordan coefficient can be written as 

(1,^1,^10,0) = -^^, (6.3) 

where / = 1 for the u quark (antiquark) and / = 2 for the d. The flavor structure 
of the T) state (as well as other isospin zero mesons) was chosen as a linear combination 
-^pA\uu) + \dd)), although in general these states are linear combinations cos((j))[\uu) + 
\dd)]/y/2 + sm(<p)\ss). The \ss) does not contribute to the 7Ti decay amplitude and therefore 
may be neglected in calculations, provided this amplitude is multiplied by a factor cos(0). 
Similarly the p (I = 1) and u, <fi (I = 0) states (J PC = 1 ) are given by 



c^ 3 Pg d 3 Pg 

(27r) 3 2£(m,p g ) (27r)32£(m, p f ) ' 

x (2vr) 3 5 3 ( Pg + p 9 - - P)-±=£ c?c? (-, -, J 3 )*J f (p„ Pq, A 9 , A 5 ) 



A,cJ 



w^" 8 ' 2' 

X iV(P) ^( m «( p <" P ^ 6 kA q / 9Cq <A ?/(?c? |0), (6.4) 

where denotes the spin-1 wave function (|5.16j) . 

The b\ (I = 1) and h\ (I = 0) mesons (J = 1 H ) have an additional orbital angular 
momentum L = 1 represented by Yjj (q) , where q is the momentum of the constituent quark 
in the meson rest frame 

C(p„ P) = A(P - 0)p„ = p, - EiShPsK + ( P '"^ (6.5) 

m q q m q q{E(m q q, -P) + m q q) 

with 

m q q = m q q(pq,P -p q ), q = q/|q|. (6.6) 

The corresponding states are given by (|6.4|) . although the spin wave function ty^g is given 
now by (ET221 . Finally, the a (I = 1) and / (I = 0) states (J PC< = 0, 1, 2++) are described 
by (|6.4|) with the spin wave function ()5.23j) . 



The states are normalized 



(P,A,/ 3 |P',A',4) = (27r) 3 2 J E(m M ,P)5 3 (P-P , )^A^/3/^ ( 6 - 7 ) 



where uim is the meson mass. That fixes the normalization constants, 

A^P) = (2^,P))-^_ ^ ^ (m ; k) ^ p _ k / j [^0(q(k,P))f 

x [Vx(m g9 -(k, P - k)/^ M )] 2 , (6.8) 



with q given in (|6.5|) and L being the orbital angular momentum of the meson. Without 
loss of generality we have taken P = Pe z . 

The orbital angular momentum wave function for a meson depends on the potential 
between a quark and an antiquark. An explicit form of such a potential is not known 
exactly and such a function must be modeled. Because of Lorentz invariance it may depend 
on momenta only through the invariant mass of a qq pair. Moreover, it must tend to zero 
for large momenta fast enough to make the amplitude convergent. The simplest choice is 
the gaussian function 

M™ q <i(p q ,P<i)/») = e-^P-P^ 2 . (6.9) 

The integrals in the decay amplitudes are not elementary and must be computed numeri- 
cally. In the nonrelativistic limit (for a large m), however, they can be expressed in terms 
of the error function. 

The free parameters of the model presented are: the quark masses m, the size parameters 
of the orbital wave functions /i and the strong coupling g. The pion decay constant /„- and 
the elastic form factor F n , defined by 

(0|^(0)|vr fc (p))=/^ fc , (6.10) 

and 

(vr''(p')|^(0)|vr fc (p)) = F n (p»+p'»)ie ijk , (6.11) 
are used to constrain the it wave function parameters (with the ir state given by (|6.1[) ). The 



axial and the vector currents are defined by 



(6.12) 



and 



^(0)=4/(0)7 M y^/, (6-13) 



with tp c f(x) given in (|7.2j) . By virtue of Lorentz invariance f n is a constant, whereas F n is 
a function of Q 2 = — (p — p') 2 - 

As mentioned previously, it is not possible to construct the wave functions with a fixed 
number of constituents in a Lorentz covariant way. Thus the current matrix elements are 
expected not to be exactly Lorentz covariant. This will be reflected, for example, in different 
values of f n obtained from spatial and time components of the axial current (rotational 
symmetry is not broken). Even if we replaced the factor E(mM,^ > ) in (|6,7|) by 1, it would 
be very difficult to find the generators of the Poincare group that satisfy the commutation 
relations. Thus, our model with the exponential orbital wave functions will not be exactly 
covariant. The resulting form factors will depend on the frame of reference. In order to 
obtain F^iQ 2 = 0) = 1, one typically employs the time component fi = and works in the 
Breit frame of reference. In this case we obtain 



\/3m f d 3 p (p° + q 



0\2 



-4 



UP) N n (P)E(m n ,P) J (2tt)3 p°q°m q g 6 **" ' {6M) 



and 



^(P,P') 



N v (P)N„(P')(E{m n , P) + E(m n , P')) J (2vr) 3 m qi m' qqP S 



d 3 p (p0 +9 0)(p0 +r 0) 



n o„o r o 



2 i 12 



x [(P • q)r° + (p ■ r)q° - (q ■ r)p° + m 2 (p° + q° + r°)]e 8 ^ , (6.15) 



where: 



q = P — p, r = P' — p, p° = E(m, p), q° = E(m, q), r° = E(m,r), 
m qq -(p, q) = [(E(m, p) + E(m, q)) 2 - (p + q) 2 ] 1/2 , 



m q g = m q g(p, q), m' qq - = m qq (p, r), 



(6.16) 



and the pion normalization constant is given in Q6.8JI . For other light unflavored mesons 
there are not enough experimental data to constrain their parameters \x. However, they are 
expected to be on the same order as fiji . 

By taking m large as compared to the //'s and Po, one obtains the nonrelativistic limit 
in which quarks are heavy. Their motion may be described by nonrelativistic quantum 
mechanics and, as we demonstrated at the end of Chapter 5, spin does not change via 
Wigner rotations. Therefore, all spin wave functions are just described by Clebsch-Gordan 
coefficients and spherical harmonics, and the spin factors in the decay amplitudes reduce 
to traces of products of Pauli matrices. All energy terms E(m, p) tend to m, whereas 
the invariant masses (|5.18|) and (|6.2U|) tend to 2m and 2m + 2?(m s ,Q), respectively. In 
the orbital wave functions, however, we must keep the next leading terms depending on 
momenta, otherwise the amplitude would become divergent: 

{Pq ~ Pq) 2 



m qq (p q ,p q ) -> 2m + 

Pq + P2 

m q qg{p q ,Pq,Q) -> 2m + m g + h E(m g ,Q). (6.17) 

In the above we have Q = — p q — Pq because the qqg state is at rest. The normalization 
constants are given in this limit by 

N 2 M (P) = 2E- 1 (m M ,P) J ^3[y L0 (q(k,P))] 2 [^(m OT -(k,P -k)//i M )] 2 , (6.18) 

where L is the orbital angular momentum of a meson. For the decay amplitudes we will 
not derive the nonrelativistic formulae from the beginning, but instead, we will go with m 
to very large values and keep only the leading terms. 



6.2 Exotic mesons as qqg bound states 

In our model a hybrid is regarded as a bound state of a quark, an antiquark and a gluon. 
Therefore, we can construct it in terms of the annihilation and creation operators. The tt\ 



state I G (J PC ) = 1 (1 + ) in its rest frame is given by 

1 f d 3 p q d 3 pg d 3 Q 



I (T \ \\ = V — [ Pq ^ 

leX[3} £X>) \cj N *t J ^2E(m, p q ) (27r) 3 2£(m, p 9 ) {2^2E{m g , Q) 

x 2(E(m, Pq ) + E(m, p ? ) + S(m fl , Q))(2vr) 3 5 3 ( P(? + p 9 - + Q) 



ex j m gqg ex' J 

x |^c 9 4' ^5 \, k\I> h)bl qXqfqC / p ^ f9C ^ XgCg \0), (6.19) 



where the spin wave function ^ qqg was given in (|5.35|) for S = 0,1,2. The orbital wave 
function ip' L depends only on the invariant mass of a quark-antiquark pair m qq and the 
invariant mass of a 3-body system, 

m qqg (p q , p q , Q) = E(rn, p q ) + E(m, p q ) + E(m g , Q). (6.20) 

Here m g is the dynamical mass of a gluon in the Coulomb gauge (arising from the strong 
interaction with virtual particles), and \ c c 9 q cq are the SU(3) Gell-Mann matrices. They 
guarantee that a hybrid meson is colorless. The gluon operators satisfy the commutation 
relations: 

[a p \c, %'x'c'] = (2vr) 3 2 J B(m 9 ,p)5 3 (p - p')5 xy 5 cd , 

[a P \c, ap'A'c'] = [Op Ac , flp'A'c'] = 0. (6.21) 

The normalization (j6,7[) leads to 

2 = 3_ (2 v-i f d 3 k dH (E{m,k)+E{m,l) + E(m g ,-k-l)) 2 (k z + l z 
ex 47T Km ex ) J (27r)3 ( 2vr )3 2E(m,k)E(m,V)E(rn g ,-k-l) (k + 1) 2 



X [^L( m qq(Kl)/Vex,rn qqg (k,l,-k-l)/ti ex ,)] 2 . (6.22) 



Hybrid mesons with other quantum numbers can be constructed in similar fashion. 

The covariant orbital wave function of the it\ may depend only on the invariant masses 



m qq and m qqg . A natural choice is a product of two gaussian functions, 

^ L (m qq (p q ,p q )/^ ex ,m qqg (p q ,p q ,Q)/^ x ) = e" ro « (p«.p*)/*& e ~ m ^ (p« >P"Q) . ( 6 .23) 
In the nonrelativistic limit the 7Ti normalization constant is given by 

^ = 4^ to^ 7 (2tt)3 ( 27 t)3 (k + 1)2 ^(m q ,(K 1)//^, m qqg (K 1, -k - 1) //W )] . 

(6.24) 



Chapter 7 



Decay of normal mesons 



The common approach to normal meson decays is based on the Po model where qq pair 
creation is described by an effective operator that creates this pair from the vacuum in 
the presence of the normal qq component of the decaying meson, as shown in Fig. 17.11 In 




Figure 7.1: The 3 Po decay of a normal meson. 

the QCD-motivated Coulomb gauge, however, the decay of a normal meson is expected 
to proceed via mixing of the qq state with the qqg hybrid component followed by gluon 
dissociation to a qq pair, as shown in Fig. 17.21 The dashed line represents the confining 
non-abelian Coulomb potential. The hybrid component of the wave function is obtained 
by integrating the qq wave function over the amplitude of the transverse gluon emission 



p 




Figure 7.2: The 3 Si decay of a normal meson via its qqg component. 

from the Coulomb line jllj . The quantum numbers P, C of such a qq state are determined 
from the corresponding conservation laws, whereas its spin may have more than one value 
(denoted in this chapter by J). 

We will be interested in estimating the size of relativistic effects in meson decays, and 
not in giving the absolute width predictions. Therefore, we can make calculations for each 
value of J separately. The relative contributions from various J and the total width may 
be obtained from the above gluon emission amplitude. Since the quark pair is emitted in 
the S = 1, L = state, this decay mechanism is also referred to as the 3 Si model. 

The Hamiltonian H of gluon dissociation is equal to 

H = J2j d 3 x^ Cl/l (x)( 57 • A c «(x))^ 2/2 (x)<5 /l/2 iAg C2 . (7.1) 

In the constituent basis used here, the single-particle quark and antiquark wave functions 
correspond to the states of massive particles with a relativistic dispersion relation, in which 
the running quark mass is approximated by a constant constituent mass m, 

f r/ 3 k 

«x) = £ J ppr^r Kk, A)6 kAc/ + „(-k, A)dl kAc/ K k ' x (7.2) 

A 

Similarly, the gluon field A Cfl is expanded in a basis of transverse polarization vectors, with 



a single-particle wave function characterizing a state with mass m 



91 

r r/ 3 k 4. 

A °' (x) = ? I iwm^xj 1 ^ A)a " + E * ( - k ' x)aX ^~ (7 - 3) 

Here g is the strong coupling constant. The Hamiltonian part contributing to the decay 
amplitude is 

" ^ J (27r)32E(m,ki) (2vr)32i?(m,k 2 ) (2ir)32E(m g , k) W d l k ~ k l + k 2 
x^Ag C2 12(k 1) A 1 )(^4(k,Aj)<-k 2 ,A 2 )6j CiAiCi/i dL k2A2C2/2 a^. (7.4) 



In this chapter we will study the decays of the p and b\ since they are dominated (almost 
100%) by a single mode, and their widths are well-known from experiment. Therefore they 
can be used to test the model presented in this work. Because we are interested in widths, 
all calculations will be done in the rest frame of a decaying meson. 

7.1 Decay p — > 2tt 

We will start from the 3 Po Hamiltonian: 

/2"C1C2 ) 

(7.5) 

with ip defined in (|7.2|) and A being a mass scale which can be fixed by the absolute decay 
width and is expected to be of the order of the average quark momentum. The amplitude 
of this mode is determined by the matrix element 

<7r(P),7r(-P)|tf|p>. (7.6) 

The corresponding states are given by Eqs. (|6.1|) and (|6.4|) . 

In this matrix element we have two operators b and two d coming from the outgoing 
meson states, whereas the decaying meson state provides one and one d) . Moreover, we 
get one and one S from the Hamiltonian. Thus each pair b, and d, d) appears twice 



and Wick's rearrangement leads to two nonzero terms. The anticommutation relations (|6.'2[) 
give Dirac delta functions that guarantee the conservation of momentum and Kronecker 
deltas acting in flavor and color space. This simplifies the integration over momenta and 
reduces summations over flavor and color indices to calculating traces of corresponding 
matrix products. 

The two terms in the above matrix element will be equal after integrating (up to a sign) 
because of symmetry, so it is enough to deal with only one and multiply the final expression 
for the amplitude by 2. A short proof of this statement will be given for the tt\ decay in the 
next chapter. Summing over color gives a factor whereas summing over flavor leads 

to a trace of a product of Pauli matrices appearing in the isospin factors ()6.3|) . 

2- 3 / 2 Tr(afa]al)e l (l0e*^)eUl^). (7.7) 

For all possible isospin channels the flavor factor is equal to ±l/\/2- For spin we obtain 



Tr 



{h + m){ p + m) (y _ m J p -j; 2m ) (/- m)(/ - m)]e\X p ) 



V^(p,r,k,l) = ^ - ;^ |1T 7 , ^— , (7.8) 

2 6 l z m q q{p, k)m qq -(r, l)m qq -{p, 1) 

where p and 1 denote respectively the momenta of a quark and an antiquark in p, and r and 
k denote respectively the momenta of a quark and an antiquark created from the vacuum, 
as shown in Fig. 17.11 We assume that all quarks are on-shell particles, i.e.: 

p° = E(m, p), r° = E(m, r), k° = E{m, k), 1° = E(rn, 1). (7.9) 

Integration over momenta gives (27r) 3 5 3 (0)A(P), where A denotes the amplitude of this 
decay, 

2A f d 3 p (E(m,p) + E(m,P - p)) 2 



ACP,\ ) = -= — / 



V6N%(P)N p (0) J (2tt) 3 E(m,p)E 2 (m,P-p) 
x W x ?(p, p - P, P - p, -p) fe(m w -(p, P - p)/^)] 2 

x M{E{m, p) + E(m, P - p))//x p ). (7.10) 



n -i 



The corresponding width is obtained from 



with Pq satisfying 



P = ^-ml (7.12) 



and Po = I Pol- The amplitude must be, according to the Wigner-Eckart theorem, of the 
form: 

A(P,X p ) = ai (P)Y 1Xp (P/P), (7.13) 
where a\ is the P-wave partial amplitude. Thus 

P " a?(P ). (7.14) 



2>2ir 2 m 2 p 



Taking P = Pe z and X p = gives 

ai(P) = y —A(Pe z ,$), (7.15) 

therefore the width is given by 

T^r e) = ^1^(^,0). (7.16) 
For large m (in the nonrelativistic limit) the trace term tends to the value 

V2V-PV(A P ), (7.17) 



and the amplitude (|7.1U|) simplifies to 

A(P,X P ) = -= ^ / T^Lfp* -P^eHXMlWr. (7.18) 

Here the normalization constants are given by ([6.18ft . and in the orbital wave functions 
we need to expand the invariant masses and energies only up to terms quadratic in 



momenta. 

Now we proceed to the 3 S± model. For the p meson the qqg component can be expanded 
in a basis of the ao, 01, «2 wave functions, all having spin 1 and one unit of the orbital angular 
momentum between the quark and the antiquark Q5.23JI . all coupled with a transverse gluon 
wave function to give the J PC = 1 state. The wave functions for the p are thus: 



]T ^ 9? (q, l qq = "Q, A„ Xg)D^*(Q)(J, X q q] 1, a\l, X p ). (7.19) 

Aqq,CT=±l 



• J.Xqq 



where tyJx qq are the ao,ai and «2 wave functions for J = 0,1,2 respectively. The 
normalized wave functions for the p are then given, similarly to those for the m <|5.33jL by: 



qqg(J=0) y 



VJ <f (q, 1 TC - = -Q, A„ A 5 -)[e*(A w -) • q][e*(Q, A,) • e(A p )], 



qqg(J=l) V 327T 



— VJ^f (q,l w - = -Q,A 9 ,A 9 -)[e*(A OT -) x q] • [e*(Q,A g ) x e(A p )], 



*qqg(J=2) 



27 ~ A . 



1607T 



^<f(q,l g9 - = -Q,A g ,A 5 ) 



xe*(Q,A fl )-[e*(Ag ? )®q]-e(A 



(7.20) 



where ^qq 9 is the spin-1 wave function (|5.16|) . Writing this function more explicitly in terms 
of the quark momenta p q and p q gives 



x VV(J)(-P</ - P$>A 9 , A p ), 



\/2m q q 



-u(p q ,Xq) 



r 



V^q-Vq 



m q g + 2m 



v(Pqi Xq) 



(7.21) 



where the gluon terms are respectively: 



^(7=0) (Qj A s , A p ) = 

^V(j=i)(Q>A g ,A p ) = 
Xe* c m (Q,X 9 )e n (X p ), 

^ti{J=2)(QAg,Xp) = 



3 
SV 
~9~ 
32^ 



K^K k 



m q q(Eqq + m g g) 

K p K k 



m q q(Eqq + m q q) 



(5 ki 5 mn )q l eT(Q,X g )e n (X p ), 

^gkm^ln ^kn ^Im^ql 
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K p K k 



m q q(E q q + m qq ) 



^km^ln _|_ gkn^lm — § k ^ ^ mn 'j 



xq l e* c m (Q,\ g )e n (\ p ), 



(7.22) 



and 



m qq = m qq (p q ,p q ), E qq = E(m qq , Q), K = -Q, 
K° = E qq + m qq , e°(Q,A g ) = 0. 



(7.23) 



As before, q = q(p g ,— Q) denotes the quark momentum in the rest frame of the qq pair 
(|6.5j) . The most general wave function will be given by a linear combination of the three 
components listed above, and the coefficients in this are provided by the qqg component 
mentioned at the beginning of this chapter. 

The Hamiltonian matrix element leads again to two equal terms. Summation over flavor 
indices gives ±l/\/2 as before, whereas for color one obtains 



1 4 

\ a \ a 

12 bc cb ~ 3" 



(7.24) 



Summation over spin gives 



(7.25) 



where 



Tr 



D 



iki 



(ft -m)(t>- m) (7" + m q ^ +2m ) (/ + m)(/ + m) 7 j 



2 3 / 2 m qq (p, k)m qq (r, l)m qq (p, 1) 



(7.26) 



and 



(7.27) 



The tensor Bj corresponds to the first term contributing to the amplitude, and the notation 
is the same as in Fig. 17.21 The functions Q7.27JI are given by: 



^(j=o)(Q' a ^ a p) = -vsr 



^(J=1)(Q' A 9' A p) 



8tt 
~9~ 
327 
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m qq (E qq + m 9 q) 



m qq (E qq + m qq ) 



- Q j Q m )e n {X p )(5 kl 5 mn ), 
q l (5 jm -Q j Q m )e n (X p ) 



n a 



x ( S km S ln _ 5 kn S lm^ 
^(J=2)(Q> A 9> A p) = 
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q k 



m q q(E qq + m qq ) 



q l (5i m -Q 3 Q m )e n {\ p ) 



^ ^km^ln _j_ gkn^lm § mn ^ 



(7.28) 



Consequently, the width can be determined from ()7.16|) . In the nonrelativistic limit 



B ij -» -yftrnti?', 



(7.29) 



and the other components are of higher order in small quantities. Therefore: 



'— mgV(A p )(^-W), 



_ m §V'(A p )(6tf + W), 
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^.m^'(A p )i(7^- W). 



(7.30) 



None of these functions vanishes. However, only two of them remain linearly independent. 



7.2 Decay b\ — ► 7tcj 

This process is a better test for this model because the ratio of the D-wave to the S-wave 
width rates is independent of the values of A and g. We begin with the b\ as a qq bound 
state and the decay Hamiltonian (|7.5|) . The amplitude of this mode is determined by the 
matrix element 

(tt(P),u;(-P)Wi>- (7-31) 

This will lead to two equal terms, as before. Summation over color and flavor gives respec- 
tively l/y/3 and ±l/\/2, whereas the spin factor is given by 



Tr 

W XuJ = — 



+ m)0fc - m)(i> - m)U- m)(>f - ^ ^ jj #-P, A.) 



2 3 / 2 m (?( j(p, k)m g g(r, l)m 9f (p, 1) 



/"■ r* 



with the same notation as in the preceding section. In the nonrelativistic limit this expres- 
sion tends to \/2(p % — P*)e**(A aJ ). The amplitude for this decay is 

A(-p \ \ ) = 2 A f d3 P (%.P)+%.P-P)) 2 

1 ' bl> y/6N K (P)N u (P)N bl (0) J (2vr)3 E(m, p)£?(m, P - p) 

x W x "(p, p - P, P - p, -p)lu 6l (p/\ P \)M(E(m, P) + E(m, P - p))/ Mbl ) 

x Vl("i 9 «?(p,P - p)//x 7r )^ L (m 9 q(p,P - p)/nu), (7.33) 

which goes in the nonrelativistic limit to 

AfP.VAfti) = -= — / - P^eHXMl^l^r 1 . (7.34) 

V ' i; v / 3mAT 7r (P)iV a; (P)Ar 6l (0)i (2tt) 3 ^ ; \ pWlvlvl \ > 

This amplitude can be expanded into the partial waves 

A(P,\ UJ ,\ bl )=Y,aL(P)Y Ll (P/\P\)(L,l;l,\ UJ \l,\ bl ), (7.35) 

L,l 

where L = or 2. The decay width is given by 

1 Pol^(Po)l 2 ^ (?36) 

with Po satisfying 

E(m p ,P )+E(m u ,-P ) = m bl . (7.37) 

This leads to 

Taking P = Pe z and X bl = X ex = gives the first equation for the two partial amplitudes, 

A(Pe 2 ,0,0) = A\\(P) = yjao(P) + \[^^{P). (7.39) 



If P = Pe±, where e± is an arbitrary unit vector perpendicular to e z , then the second 
equation for the partial amplitudes is 

A(Pe ± ,0,0) = A ± (P) = J]-ao(P) - \[^a 2 {P). (7.40) 

V 47T V 87T 

Finally, one obtains: 

r fcr e) = 72^r^ll( P o)+2^(P )] 2 , 

<^ e) = ^iMPo)-A ± (P )?. (7.41) 

Now we move to the b\ treated as a gluonic bound state. The qqg wave function with 
J PC = 1 H , 1=1 quantum numbers requires the qq to have the 7r or 7T2 quantum numbers. 
The corresponding, total wave functions are given by 

Klkj)^ = E W = "Q. V \-)4 1 s t*(Q)(^ 1, <r|l, A fel ), 

A<jq,cr=±l 

(7.42) 

with ty'g q9 being the ir (7^) gc/ wave function for J = (J = 2). The normalized spin wave 
function is thus given by 



^^(J) =E^-(q'-Q' A 9 ' A 5X(j)(Q,A,A 9 ,A /C) ), (7.43) 

A 



where 



C(J=o) = vJ[e*(Q,A fl )-e(A 6l )] 



C(j=2) = • [e*(Q, A 9 ) ® e(A 6l )] • q, (7.44) 

with Q being the gluon momentum and q = q(p 9 , — Q) being the relative momentum in 
the qq pair (|6.5|) . The spin factor is given by (|7.25|) , but with a different tensor : 



Tr 

B« = — 



(/s- m)(i> - m)(/- m)(V - m J(r,i)+2m ) (/ + m )7 J 



2 3 / 2 m^(p, k)m Q g(r, l)m g5 (p, 1) 



(7.45) 



The notation used above is the same as in Fig. 17.21 In the nonrelativistic limit, however, this 
tensor tends to the same values as the B w defined in l|7.26|) . The widths can be calculated, 
as in the 3 Po model, from (|7.41|) . 

7.3 Numerical results 

In a simple constituent quark model one assumes that the mass difference between the 
mesons tt and p arises only from spin. Therefore, we can write 

rriM = mM + k(si ■ s 2 ), (7.46) 

where M denotes either meson and fhu is its "averaged" mass. Using the identity s(s + l) = 
+ 1) + S2(«2 + 1) + 2si • S2, we can solve for %. We have s = for the n, and s = 1 
for the p. We have also s% = S2 = 1/2. Substituting =140 MeV and m p =770 MeV 
gives ffiM =612 MeV. Thus m u = rrid = fhM/2 =306 MeV. A similar relation can be used 
for the K and K* (decays to strange mesons will be described in the next chapter), leading 
to m K =792 MeV and m s = m K - m u =486 MeV. 

The averaged mass of the tt and p mesons should rather be used instead of their physical 
masses in the normalization constant (|6,8|) . In similar fashion, the averaged mass of the K 
and K* should be used in ()8.45|) . This procedure can also be followed for the b\ and aj 
mesons, or for the h\ and fj mesons. In this case, however, there is an additional term 
proportional to the spin-orbit interaction, S ■ L. The averaged masses for the b\ and fx 
are found to be close to their physical values. Therefore the physical values will be used in 
normalization for these mesons. 

The weak decay constants (|6.14j) and (|8.46|) can be used to fit the parameters p^ and 
Pk- Because our model is not exactly Lorentz covariant, the weak decay constants become 
functions of the meson momentum and we choose them to be equal to their experimental 
values at rest. Thus, setting f n (0) = 93 MeV and f K (0) = 113 MeV leads to p n = 221 MeV 
and pk = 275 MeV. The momentum dependence of in our model for m = 306 MeV and 
p^ = 221 MeV is presented in Fig. 17.31 which shows the difference between the rest-value 
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Figure 7.3: The pion weak decay constant f n as a function of the pion momentum p. 

and the infinity- value at the level of 20%. 

The strong coupling constant at this scale is approximately g 2 = 10, and for the effective 
mass of the gluon we will take m g = 500 MeV, following what we said in Chapter 3. In 
Fig. 17.41 we present F 2 (Q 2 ) calculated with the same wave function parameters and com- 
pared with data [Sj. The agreement is good for small momentum transfer, whereas the 
discrepancy for larger Q 2 indicates a missing, high momentum component of the wave func- 
tion. Actually, we could constrain the quark mass m from and F n but it appears that 
both quantities are not too sensitive to m. Thus its value taken from the averaged mass of 
the 7T and p mesons works pretty well. 

For the p and b\ decays we will choose the values of all parameters p to be equal to p n . 
The numerical predictions for the p widths with A = p n are presented in Table [7T1 As we 
discussed at the beginning of this chapter, these numbers correspond to each qq spin value 
in the hybrid component of a decaying meson. The experimental value of the width for 
p — > 2ir is 149 MeV |78j . This number can be used to fit the free parameter A, the coupling 
constant g, or the size parameter p! ex which need not to be of the same order as In 
order to do so, however, we need to know how each J contributes to the total p spin wave 
function. 

The numerical predictions for the b\ widths with A = p n are presented in Table 17.21 
The experimental value of the total width for the process b\ — * txuj is 142 MeV, and for 
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Figure 7.4: The pion electromagnetic form factor F% as a function of the momentum 
transfer square Q 2 . 
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Table 7.1: Relativistic (nonrelativistic) widths in MeV of the decay p(770) — ► 2ir for the 
3 Po model, and for the 3 Si model with p = aj + g and J = 0, 1, 2. 

the ratio of the D-wave and S-wave width rates is 0.08 |78j . Our predictions give a value 
less than 0.02 for this ratio in the 3 Si model, and close to 4 for the 3 Pq decay. Therefore 
the real mechanism should lie somewhere in between, although the 3 Si mechanism gives 
more accurate result. However, in the 3 Po model the D/S ratio is very sensitive to the free 
parameters p and for p =400 MeV one obtains this ratio on the order of the experimental 
value ^S]. This value of p, however, is inconsistent with the weak decay constant and the 
elastic form factor for the pion. 

The qqg wave function component of the b\ wave functions used here is that of Eq. I|7.42|) . 
corresponding to a qq pair with the tt quantum numbers. For a qq with the TT2 quantum 
numbers, the numeric value for the width is much smaller than 1 MeV for the S-wave 
and approximately 1 MeV for the D-wave. The ratio D/S is respectively 230. In the 
nonrelativistic limit we obtain similar results. 

We observe that treating the b\ as the tt2 + g state increases dramatically the D/S 
ratio. Therefore this may be an important component of the wave function. Relatively 
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Table 7.2: Relativistic (nonrelativistic) widths in MeV of the decay &i(1235) — > ^6^(782) 
for the 3 Po model, and for the 3 Si model with b\ = tt + g. 

small values of the decay widths of b\ with 7T2 quantum numbers (L=2) compared to those 
of b\ with pion quantum numbers resemble the situation for the process t\\ — > 7r&i, whose 
D-wave width was small compared to that in the S-wave. 



Chapter 8 



Decays of tt\ 



In this chapter we will study the main subject of the presented work, i.e., a completely 
relativistic decay of the exotic meson tt\. In experiment we observe that most hadronic 
decays involve a minimal number of final state particles. This requires a small number of 
transitions at the quark level. Thus, an exotic meson is expected to decay into two normal 
mesons. In a constituent quark model the transverse gluon in the tti dissociates into a 
quark and an antiquark, and the two resulting quark-antiquark pairs rearrange themselves 
into two mesons, as shown in Fig. 18.11 The possible decay modes of the tti for m 7Tl = 1600 
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Figure 8.1: The 3 Si decay of a hybrid into two mesons. 
MeV are listed in Table I87T1 where L is the angular momentum between outgoing mesons. 
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Table 8.1: Possible decay modes and the angular momentum between the outgoing mesons 
L for the iri with mass 2.0 GeV. 

For brevity, we did not put there the decays into the antiparticles of the corresponding 
strange mesons. Hereinafter, such decays are understood to be included and have the same 
widths. In the following sections we will discuss the most important modes and calculate 
the corresponding widths. 



8.1 Decay of m into tttj and nbi 

The amplitudes for these modes are related to the matrix element 

(TrfP^MC-PJIfllTTi), (8.1) 

where M is either rj or b\. From rearranging the annihilation and creation operators one 
obtains two nonzero terms contributing to the total amplitude. We will show below that 
they are equal. Summing over color gives a factor 4/3, as for decays of the qqg component 
of a normal meson. Summing over flavor leads to a trace of a product of Pauli matrices 
appearing in the isospin term (|6.3|) . 

2- 3 / 2 Tr(a?al)e i (II*)e*(lZ) (8.2) 



TO 



or 

2-WTr(ora?o%)e i (I?y j (IZy k (#), (8.3) 

respectively (for all allowed isospin channels these factors are again equal to ±1/V2). The 
spin factor is given by (|7.25j) with S = J and X p replaced by \ ex , where the tensor B M was 
introduced in (|7.26|) and ififi* are given by Q7.27JI . The notation is shown in Fig. 18.11 Using 
(j5.36j) we can write: 

Vvf (5=0) = -{^9%V k - Q j Q k )Q l e l (X ex ), 

* \8TTl^ m q q(E qq + m qq ) i 

^ (5=2) = ^W (5=1) -^ (5=0) ). (8.4) 
In the terms above we used the following notation: 

mqq = m qq (p, -P-r), Q = P - p + r, E qq = E(m qq ,Q), K = -Q, K° = E qq + m gq , 

(8.5) 

and assumed all quarks being on-shell particles. 

Integration over all momenta gives for the first term (27r) 3 <5 3 (0)Ai, where the amplitude 
A\ is given by 

AS) cp \ \ n _ M 3 P d3r (E(m,p) + E(m,k))(E(m,r) + E(m,l)) 
1{M) ( ,*M,*ex) J (2tt) 3 (2tt) 3 4E(m,p)E(m,-k)E(m,r)E(m,l)E(m g ,Q) 

x (E(m,p) +E(m,l) + E(m g , Q))ip' L (m qq (p, V)/fJ, ex , m qqg (p, 1, Q)/// e J 
x[N^P)N M (P)N ex ]- l ^ L {m qq {p^)/^ L {m q ^ 

Here, M denotes the second meson (rj or b\), J is its total spin, q = q(r, — P) is given in (|6.5[) 
and P = |P|. If M = r), then the above expression must be proportional to e l {\ ex )P l (or 
Y\\ ex {P / P)), being a vector function of the vector P (the outgoing mesons in the P-wave). 



rr A 



However, if M = b\ then the integral (|8.6|) is proportional to 



e l (\ ex y*(\ M )[P l P J + f(P)5% (8.7) 

and can be represented as a superposition of spherical harmonics corresponding respectively 
to I = and I = 2 (S-wave and D-wave). 

(S) 

The second term in the Hamiltonian matrix element A 2 is obtained from the first one 



by interchanging 



p «-► r, P <-> -P, iP 4 2) 



everywhere in ()8.6|) . including k, 1, Q, except in the third line. For tv\ — * irr] we can write 
schematically: 



Ai(p; 

x (P - p + r 
A 2 (P 



y d 3 p ( i 3 r/ r) (mgg(p,P - p))/ 7r (m gf (r,-P - r))/ ea .(m g g(p, -P - r)) 



y d 3 pd 3 rf v (m g g(j),'P - p))U(m q9 (r, -P - r))/ M (m ?f (r, P - p)) 
x(-P-r + p). (8.9) 

We have not included here the energy and trace factors because they do not change under 
the above symmetry. Using m qq -(p, r) = m q q(r, p) = m q q{— p, — r) one can show that both 
terms are equal (for the same A m and X ex ), 



A 2 (P) = -Ai(-P)=Ai(P). (8.10) 

Thus, we have A = A± + A2 = 2A\. In similar fashion one can prove this equality (up to a 
sign) for all decays of any qqg particle into two mesons. 

If [A„ = fi n , then each term (Ai or A2) is a product of a part that is symmetric under 
interchanging p *-* r + P and a (vector) part that is antisymmetric. In this case the hybrid 
will not decay into ir and 77. Neither can it decay into two pions because of a minus sign 
from interchanging the Pauli matrices in the isospin factor that makes both terms cancel. 
There is the same minus sign for the 7Ti — * nbi, but now the amplitude is a scalar (L = 0) or 



a tensor (L = 2) function of P and again Ai(P) = ^(P). However, in this case p bl = Pn 
does not imply A = because the orbital wave functions of these mesons are different. 
Thus we find that the 1 ^ isovector does not decay into identical pseudoscalars. This is a 
relativistic generalization of a symmetry found in other nonrelativistic decay models |54j . 
One might think that a similar symmetry could occur for the decay p — > 2tt. In this case, 
however, the amplitude does not vanish because its vector part comes from a spherical 
harmonic associated not with the gluon momentum, but with the relative momentum in a 
qq pair. This harmonic is asymmetric under interchanging p *-* r + P and the amplitude 
for the p — ► 2-7T does not vanish. 
The width is given by 

ff _ i fli^QE 

32tt z mj x 

with Po satisfying 

E(m n ,P ) + E(m M ,-P ) = m ex (8.12) 

and Po = | I*o I - The amplitudes must be, according to the Wigner-Eckart theorem, of the 
form: 

A^(P,X ex ) = a^(P)Y 1Xex (P/P), 

A^(P,X bl ,X ex )=^al h (P)Y u (P/P)(L,l,l,\ bl \l,X ex ), (8.13) 

L,l 

which leads to Eq. (|7.38|) with m bl replaced by m ex . The amplitudes A(P, Xm, X ex ) given 
in ()8.6j) are multiplied by 2, since there are two equal terms. The expression for the tt\ — > 
ttt] width is similar to (|7.16|) . whereas for the tti — > ixb\ we can use the "parallel" and 
"perpendicular" amplitudes introduced for the decay &i — ► ttlo. The final results are: 

Up-wave) _ ^0 A 2 I p r,\ 

L ™> ~ 24vrmL ° " J ' 

v (s-wave) _ Pq f /t 7rb i f p \ . 9 ATtbi / p \i2 



ex 



e.e 



The 7r and r\ mesons have the same quantum numbers (except isospin) so fj, n and \x n 
should be almost equal. This equality is not exact because the SU(3)f is only an approxi- 
mate symmetry and there is a contribution of the ss in r]. Therefore the amplitude for the 
7Ti — > irq should be close to zero, and of two channels irr], nb\ the latter will be favored. 
However, the free parameters \x need not to be close to each other for decays of the tt\ 
into two mesons with different radial quantum number, which would make such channels 
significant. 



8.2 Decay of m into irp, irfi, 71/2, rjai and 7702 

For these channels the procedure is analogous to that in the preceding section. The ampli- 
tude of this process is given by the matrix element (Mi(P), M2(—~P)\H\iri) which is again 
a sum of two terms, and Mi and M2 denote the two outgoing mesons. The color and flavor 
factors are again 4/3 and ±l/\/2, but the spin factor is given now by 



j ' 



(8.15) 



where (for the first term) 



C M = _[ 2 3/2 mM _(p jk ) m ^( r) i) mM _( P) i)]-irr[(^ + m)(y 



m q g(p, 1) + 2m 



V-m)(f 



m qq -(r, 1) + 2m 



) (/ + m) r f( / k — m)7 £ 



(8.16) 



and 



^' (5) (A, A e:c ) = ]T (Q, Xg, \ ex K(Q, W* (-P, A) 



(8.17) 



or equivalently: 



V J(S=0) 
V J(S=1) 



V J(5=2) = 



^g*(6 jk - Q j Q k )Q l e l (\ ex y*(-P,\), 



9l 



Kfj,K k 



m q q(E qi j + m qq ) 



(P l -Q^Q l )Q K e l (X ex )e^(-P,X), 



vn l \n k J 



3 - ■ vj(s=i) _ 2 , „j(s=oh 
3^ J ' 



(8.18) 



The notation used above is the same as in (|8.5[) and Fig. 18.11 
The amplitude for the first term is given by 



A[f M) (P,\ M ,\e X )=Y, [ 
X,l J 



d 3 p d 3 r (E(m,p) + E(m,k))(E(m,r) + E(m,l)) 
(2?r) 3 (2vr) 3 AE(m, p)E(m, k)E(m, r)E(m, l)E(m g , Q) 



x(E(m,p) + E(m,\) + E(m g , Q))ip' L (m qq (p, \)/p, ex , m qq - g (p, 1, Q)/pf ex ) 

x [N n N M N ex ]- l ^ L (m qq {p, k)/p n )^ L {m qq (r, l)/// M )Y L * 9 _ (M)J (q/|q|) 

x (1, A; L qq{M) J\J M , X^gl^C^fyiX, X ex ), (8.19) 

where M denotes the second meson [p with L qq = or fj with L qq = 1), J is its total spin, 
q = q(r, — P) was defined in Q6.5|) and J = 1, 2. Again one can show (for the same A's) 



4 S) (P) =4 5) (P)- (8-20) 

If //p = /Utt, then the amplitude of the decay into np is not symmetric under interchanging 
p <-> r and does not vanish, unlike for the nr] case. Therefore this mode may be significant. 
The same holds for the irfj channels. 

From the Wigner-Eckart theorem we have: 



A^(P, X p , X ex ) = al P (P)Y Ll (P/P)(L, I; 1, A p |l, X ex ), 

L,l 

A^-'{P,Xf J Ae X )=Y. a l fj ( P ) Y Ll^lP){L,hJAf J \lAe X )- (8.21) 
L,l 

In formulae (|8.21j) and ()7.38|) one must use the ^4(P, A's) given in (|8.19() . multiplied by 2. 
The mesons in the irp channel go out in the P-wave, for the nfi it is either the S-wave 
or D-wave, whereas in the irf2 they go out only in the D-wave. Taking P = Pe z and 
A p = X ex = +1 leads to 

a7(P) = \f^ A*r(Pe z ,l,l), (8.22) 



and therefore 



rfr wave) = Y^r A U p ^ L !)■ ( 8 - 23 ) 



In the above Pq satisfies ()8.12|) with tum = m p- For decays into nfj one may follow 



the procedure with the parallel and perpendicular amplitudes described in the preceding 
section, or just use the general formula that is equivalent to (|8.21|) . 




and substitute it into (|7.38f) with replaced by m ex . Here, Po satisfies (|8.12(l with 
rriM = m fj- All results for the decays iri —* irfj are valid also for the rja j modes (with 
different normalization constants, parameters [i and Pq). In this case, n must be replaced 
with 7] and / with a. 



8.3 Decay tti —> puj 

The amplitude of this process is given by the matrix element (p(P), w(— P)|iJ|7Ti), being 
again a sum of two terms. The color and flavor factors are the same as before, and the spin 
factor is equal to 

D^Jllj, (8.25) 

where (for the first term) 



D^Pi = [2 3 / 2 m gg -(p,k)m (?9 -(r g2 ,l)m (?9 -(p,l)]- 1 Tr[(/ : -m)(7 i - 



x {p + m)(^j 



m qq (p, 1) + 2m 



(/-m)( 



r 



m qq (r, 1) + 2m 



m qq (p,k) + 2m' 
(/ + m)7 i l, (8.26) 



and 



E VV(S)(Q> V A e *H(Q, A 9 )e-(P, A p )e"*(-P, A, 

A 



or equivalently: 



.27) 



^pKS=o) = _ I gk^k _ QjQk)Ql e l (\ ex y*(P, A p )e"*(-P,A W ), 



f P j(S=l) = _, ±_ 



9u 



m qq (E qq + m qq ) 



(5 jl - Q j Q l )Q k e l (X e 



xe-(P,A p )e^(-P,A w ), 

vpj{S=2) . 



.28) 



The notation used above is the same as in the preceding sections. 

Again one can demonstrate that the second term in the amplitude is equal to the first 
one. Therefore the total amplitude is 



A^ s \P,X p ,X L ,,X ex ) = 2 J 



d 3 p d 3 r (E(m,p) + E(m,k))(E(m,r) + E(m,l)) 
(2vr) 3 (2tt) 3 AE(m, p)E(m, k)E(m, r)E(m, l)E(m g , Q) 
x (E(m,p) +E(m,l) +E(m g ,Q)) ■ i/j' L (m qq -(p, l)//J, ex , m qqg (p, 1, Q)/// e J 

4 1 
'372 



x [NpN^N^ijUm^p, k)/^)feK-(r, \)/^)g -—D^ p ^ p r ( 8 - 29 ) 



If \x p = fj,u, then the symmetry of the orbital wave functions causes A = and the hybrid 
will not decay into p and u>. Because both parameters \i are expected to be on the same 
order, the pco mode will not be favored. The width is given by (j7.38j) with replaced by 
m ex , where 



'(P)= ]T f dn P A^(P, X p , X^, X^Y^P/P)^, X p ;l, XJJ' A')(L,l; J' , X'\h X ex ), 
L,l,J>,\> 

(8.30) 

and Po is obtained from E(m p , Po) + E(rn w , — Pq) = m ex . In this process we have either 
L = 1 or L = 3. 

8.4 Decay into strange mesons 

The above results can be straightforwardly generalized to the case where m q and nig are 
different, for example to decays into mesons with one strange quark (/ = 1/2). The spin 
wave function for a quark-antiquark pair in a J p = 0~ state is 



WqqQ-q, k, <V Xg) = —f=z. T. \ M m q , \ q , X q )"f 5 v(m q , \ q , Xg), (8.31) 

V2m qq (l q ,l q ) 

and the K states (understood to be both K and K) are given by 



l , ^ , . . , l 



all A,c 



(2n) 3 2E(m q ,p q ) (2ir) 3 2E(mg,pg) 

—^ScqCq^qqiPqi Pq, \i Xg) Jyfp\ 



X (27r) 3 (5 3 (Pq + Pg ~ P)^=5 CqCq fygg(p q , P q , Xg, Xg ) — ^— - l\) L (jXl qq ( Pg , Pg ) /jj) 



X 6 p,A,/,c 9 4> f A ? / fC? l )> 



.32) 



where the operators satisfy the anticommutation relations: 



{ ft pA/c> fe p'A'/'c'} = (2TT) 3 2E(m q ,p)5 3 (p - p')5x\i5 fr 5 cc ,, 



VA7 

{dp\fc,d p >\>f> c r} = (2vr) 3 2 J E(m 5 ,p)5 3 (p - p')5x\>6ff6 c 



.33) 



The invariant mass is defined as 



m q q = m q q(l q , Ig) = J (E(m q , \ q ) + E(m q , l q )) 2 - (\ q + l q ) 2 



S.34) 



and the "modified" invariant mass is 



rh qq = \Jm 2 qq - (m q -m q f 



.35) 



The flavor indices f q and f q give four combinations: us, ds, su and sd, corresponding to 
appropriate mesons. 

If in (jST52Jl *$> qq is given by 



-u(m q ,\ q ,\ q ) 



m _ 

L q l q 



m qq + m q + m q 



V(m q , lg,Xq)€ fl (lqq,\qq), 

(8.36) 



then one obtains the K* states (J p = 1 ). For a qq pair with spin S qq = and the orbital 
angular momentum L = 1 we have 



^ qq V-qi ^-qi ^qi ^q) 



V2m q q(lq,\q) 



(m q , l q , X q )^v{m q , l q , X q )Y 1Xq .(q), (8.37) 



which corresponds to a K\ meson (J p = 1 + ). For S qq = 1 we get 



1 



¥ f (l q ,l q ,X q ,Xq) = u(m g , lq, A g ) 

x e M (l g5 , A)y u (q)(l, A; l,l\J, X qq ), 



r 



lq lq 



m qq + m q + m q 



v(m q , lq, Xg) 

(8.38) 



o -i 



where q was defined in (|5.24[) . The last wave function characterizes mesons Kq, K\ and 
K* (J p = 0, 1, 2+). In a similar manner one can obtain states with higher orbital angular 
momenta L. 

The decay of tt\ into two strange mesons proceeds through the dissociation of the gluon 
into a pair ss. For the processes tt\ — > KK* and ix\ — * KK\ the Hamiltonian matrix 
element is H = (2tt) 3 5 3 (0)A (now there is only one term). The decay amplitude is given by 



E 



d 3 p d 3 r (E(m,p) + E(m' ,k))(E(m' ,r) + E(m,l)) 
(2vr) 3 (2tt) 3 4E(m, p)E(m', k)E(m' ,r)E(m, l)E(m g , Q) 



x (E(m,p) + E(m,l) + E(m g , Q))if>' L (m qq (p, I)/ Hex, m qqg (p, 1, Q)/fi' ex 

x [NKNuNe^-^Limqqip, k)/fi w )tp L (m qq (r, l)/^M)^L* M «(q/|q|) 

4 1 

x (1, A; L M , 1\Jm, ^M)g^—^(spin){X, \ ex ). 



(8.39) 



Masses of quarks u,d (assumed to be equal) and s are denoted respectively by m and m', 
whereas M stands for either K* or K±, and S = S q q g . The spin factor is given by 



AO 



?.40) 



for kaons with S qq = 1, and 



5.41) 



if Sgo = 0, where 



Tr 



2 3 / 2 m 9 ,j(p, k)m g q(r, l)m (?5 (p, 1) 



.42) 



and 



C"^' = -[2 3 / 2 m w -(p,k)m M -(r,l)m w -(p,l)]- 1 Tr 



■ + m) 7' 



m g g(p, 1) + 2m 



(/-"») (7" f~T\ : . 7 )(/ + m')7 J (^-^')7 5 

V m qq (r, l) + m + m ' 



5.43) 



The notation used above is the same as that in (|8.5|) and Fig. 18.11 with the invariant mass 



Q8.34JI and some modifications due to difference between the masses of quarks: 



p° = E(m, p), r° = E(m, r), A; = E(rri , k), 1° = E(m, 1). (8.44) 



The normalization constants for strange mesons are 

Ar9 ,„ N , i /" d 3 k (£(m,k) + E(m'P - k)) 2 . . „ 



x [fe(m, ? (k, P - k)/ w )]', (8.45) 

where q = q(k, P) is in (|(i.5|) . The kaon weak decay constant can be used to constrain the 
value of hk- 



JK{ 1 N K (P)E(m K ,P) J (2tt)3 p°g°m 9f (p,P-p) • ^ J 



8.5 Nonrelativistic limit 

For the decays it\ — > 7T7/, 7r&i the behaviour of the tensor (|7.26|) is given by Eq. (|7.29|) . 
Therefore the spin factor B^ipj^ is: 



5 iC ^v^ m ^ e ( A -)' 

^^f 2 ^--^^f°). (8.47) 



From the above it follows 

r (s= i)^0, (8.48) 
1 (5=0) 1-3 

The amplitudes ^^(P) and yl T 6 1 (P) given by 1)8.13(1 are in this limit (for 5 = 0) 

A„ M (P) = -2--^- ^ NAp)NM{p)Nex J (2^)3 (^)3^o(q/|q|) 

ex i m qqg (p,l,Q)A4), (8.49) 



where q = p q2 + P/2, X ex = and k, 1, Q were defined in (|8.5|) . As before, M denotes 
r) or b\, with J = or 1. In the nonrelativistic limit only the 5 = and 5 = 2 hybrid 
components of the ~k\ state contribute to its decays into irr] and irbi. Therefore one may 
expect the 5 = 1 component not to be too large for a fully relativistic case. 
For 7Ti — > 7T/3, 7r/i : 2 we have 

C**J -> ^ime oikj (8.50) 
and the other components are much smaller. This results in 



C^%f - 0, 



Therefore we have 



\ S =o) 0, tT" 

1 (S=l) 



r (5=2) 9 



13 



.52) 



The amplitudes A np (P) and ^4.^ 2 (P) given by (|8.21l) are in this limit (for 5 = 1) 



,™ o 4 1 Hi~2m + m„ g f d 3 p d 3 r - . .. .. 

^ M(P)=2 37f m V4^ m 2 mg N n (P)N M (P)N ex J (2^ (2^ Q ^«-) o(q/|q|) 

(p,l,Q)//4) 

x (CG) M) (8.53) 

with the notations from Q8.5|) and M standing for p {L qq - = and (CG) = 1) or f\ 2 (L g g = 1 
and (CG) = l/y/2). In the nonrelativistic limit only the 5 = 1 and 5 = 2 qqg component 
of the 7Q state contribute to its decays into np and 7r/x,2- Therefore one may expect the 
5 = component not to be too large for a relativistic case. Similar conlusions can be made 
for the 7701,2 modes. 

Finally, for ix\ — * puo we have 

D m ^ ^2 m {5 ij 5 kl - 5 ik 5 jl + (8.54) 



Ci A 



(the other components are much smaller), which gives: 

^3~ 



D ^ P Al~ P P -» V ^{Q^ - g J ^)e"(Ap)e J *(A w )e fc (A ea; ). (8.55) 

The amplitude A puJ {P) in this limit is described by formula 1)8.49)1 (for 5 = 0) or 1)8.53)) (for 
5 = 1), where tt is replaced by p and M = to. 

We can straightforwardly understand the difference in amplitudes coming from the spin 
wave function. If we assume = m p , p v = p p and = mf 1 = m/ 2 , pb± = pj L = pf 2 
(the second condition for masses is satisfied to a good approximation), then comparing 
with (IH331 and (EH!) with (EH3)l gives 

4 — - a _> r — lr 

Ar/l = ^7/^ 7rb l ~* = g 1 **!' (8.56) 

Her©, Afpny Ajrfj-^ were taken for 5 = and A™, 2 for 5 = 1. The relation between 
A vr] and A^ (or between and A 1T f 1 2 ) is more complicated and depends on the orbital 
angular momentum wave functions ipL and ip' L . If p p = p n , then in the nonrelativistic limit 
7i"i would not decay into irp since the corresponding amplitude becomes proportional to that 
of the irrj mode. Therefore the width of the irp mode is expected to be much smaller than 
that of 7r&i, assuming the parameters p r ho and p n are very close to one another. 

Analogous calculations can be made for decays of 7Ti into strange mesons. The ampli- 
tudes A K k i (S q q = 0,1) behave similarly to A^ and A^fr, whereas A K k* behaves like 
A np . Therefore the former will be dominant and the latter is expected to be much smaller. 
If px* = Pk> then in the nonrelativistic limit A K ^* = 0. 



8.6 Numerical results 

In Tables 18.21 and 18.31 we present the widths for various decay modes containing radially 
ground-state mesons. The numbers in parentheses correspond to calculations using the 
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Table 8.2: Relativistic (nonrelativistic) widths in MeV for the tt\ decay modes with mass 
1.6 GeV, for all possible values of the total intrinsic spin S of tt\. 
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Table 8.3: Relativistic (nonrelativistic) widths in MeV for the tt\ decay modes with mass 
2.0 GeV, for all possible values of the total intrinsic spin S of tt\. 

nonrelativistic formulae, and S denotes the total spin of the qqg component in the ir\ wave 
function. As before, we are interested in relativistic effects in the tt\ decays. Thus it is 
sufficient to calculate the amplitudes for each value of S separately. The procedure for 
obtaining the relative contributions for each S was given at the beginning of Chapter 7. 
The values of parameters p for all unflavored mesons and n± were taken equal to fj, n , and 
for all strange mesons equal to fix- This makes the width values for the modes irq, nrf' and 
puj identically equal to zero. 

The mode np was expected to be rather small because in the nonrelativistic limit the 
corresponding amplitude was proportional to that of the mode 7rr/, assuming fi p = fi v . We 
also assumed \i p = p w . Combining both equalities led to this conclusion. We also found 
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Figure 8.2: Relativistic, nonrelativistic, and semirelativistic widths for the S-wave decay 
it i — > 7r6i as a function of the quark mass m, for a 7Ti with mass 1.6 GeV. 

that the qqg component of tt\ state with S = did not contribute to the hybrid decay in this 
limit identically. For the relativistic case, however, we observe that the largest contribution 
in the np channel comes from the S = component. Our other predictions are confirmed 
by numerical results, i.e., the S = 1 component for the nbi mode and the S = component 
for the 7r/i mode are negligible in the relativistic model as well. 

The decay mode irbi is a dominant one. In the following, we will analyze how its widths 
depend on the free parameters, except the coupling constant g which just multiplies the 
amplitude and \x ex whose equality to fj, n seems to be justified. Because the D-wave value in 
this case is small as compared to the S-wave, we will only deal with the latter. Furthermore, 
in the following plots we will assume S = since S = 1 gives negligible contribution (in the 
nonrelativistic limit it is zero) and thus the widths for S = 2 are proportional to those of 
S = (|8.48|) . The default values of the free parameters are those from Table 18.21 

In Fig. 18.21 we compare relativistic and nonrelativistic predictions for the width of the 
decay %\ — * ixb\ as a function of the mass of the light quark m. It also shows the semirel- 
ativistic values which include relativistic phase space and orbital wave functions, but with 
no Wigner rotations. For larger values of m all three curves R, NR and SR converge, as it 
should be. The corresponding ratios NR/R and SR/R are shown in Fig. 18.31 

Let us proceed to analyze the contribution of the other free parameters. We show only 
relativistic and semirelativistic rates because the nonrelativistic ones correspond to different 
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Figure 8.3: Width ratios NR/R and SR/R for tt\ — ► irbi, as functions of the quark mass 
m, for a tti with mass 1.6 GeV. 
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Figure 8.4: Decay widths for tti — > irbi as functions of tt\. 



orbital wave functions and we are interested in corrections arising from Wigner rotations. 
The dependence of the width of tt\ — ► irbi on the mass of tti is presented in Fig. 18.41 We 
observe that Wigner rotations contribute more for larger values of m ex , and at 2000 MeV 
the ratio SR/R is already about 50%. We also see that the maximal width occurs for m ex 
in the range 1500 — 1600 MeV, i.e., for the values close to those experimentally reported. 
Thus, our predictions for the width are an upper limit with respect to m ex . 
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Figure 8.5: Decay widths for tt\ — ► nbi as functions of the parameter iz^, for a 7Tj with 
mass 1.6 GeV. 
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Figure 8.6: Ratio SR/R for wi — > Trbi as a function of the parameter /i^, for a ix\ with 
mass 1.6 GeV. 

In Fig. 18. 51 we show the width dependence on fi^. If the value of //^ is equal to /%, then 
the width is approximately at its maximum. Probably it is a coincidence, but it tells us 
that our prediction is an upper limit for this width with respect to \i\, x . The corresponding 
ratio SR/R is given in Fig. 18.61 At masses around 1 GeV it already behaves like a constant. 

There is one parameter whose choice was not justified, [i ex i. This quantity divides the 
invariant mass of the qqg system, and there is no reason why it should be equal to fj, n , as 
we assumed. Its value may be obtained from decays of normal mesons such as p or b\. In 
Fig. 18.71 we show how the width of 7ri — ► irbi as a function of [i ex ' ■ The ratio SR/R is given 



in Fig. 18.81 and we see that it is close to a constant, i.e., the choice of fi ex i is crucial for the 
widths but not for relative relativistic corrections. 
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Figure 8.7: Decay widths for 7Ti — > irbi as functions of the parameter fj, ex ', for a 7Ti with 
mass 1.6 GeV. 
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Figure 8.8: Ratio SR/R for %\ — ► i:bi as a function of the parameter fi ex >, for a i\\ with 
mass 1.6 GeV. 

Finally, in Fig. I8.9l we show the width dependence on the effective mass of the gluon m g , 
and in Fig. !8.1Ul the corresponding ratio SR/R. The latter displays clearly that the ratio is 
essentially independent of m g . We observe that the relative relativistic corrections arising 
from Wigner rotations are insensitive to the free parameters referring to the tti at rest, i.e., 
m ex i and m g . 
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Figure 8.9: Decay widths for 7Ti — ► irbi as functions of the effective mass of the gluon m g , 
for a tt\ with mass 1.6 GeV. 




0.5 ' 1 1 1 1 1 1 1 1 1 1 

200 400 600 800 1000 

m, [MeV] 

Figure 8.10: Ratio SR/R for 7ri — * irbi as a function of the effective mass of the gluon m g , 
for a tti with mass 1.6 GeV. 

Now we move to the decays into n and various rj mesons. Since the amplitude vanishes 
for jj,^ = we want to observe what happens if fj,„ ^ fi w . For brevity, we will consider only 
the relativistic case. In Fig. 18. ill we show how the width of the decay tv\ — > irr] depends on 
the value of /U„, while the other parameters are the defaults. The biggest value approached 
is only a few MeV. That confirms our previous predictions and is in a good agreement 
with other models [31]. The same dependence for the decay tt\ — * 7rr/(958) is displayed in 
Fig. I8.12( and the maximum is around 10 MeV. 
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Figure 8.11: Decay width for 7ri — > nr] as a function of the parameter fx„, for a tt\ with 
mass 1.6 GeV. 




Figure 8.12: Decay width for tti — > 7r?/ as a function of the parameter for a 7Ti with 
mass 1.6 GeV. 

We see that the maximum of the width occurs when fi^ is approximately equal to 100 
MeV. In Fig. 18.131 and Fig. 18.141 we show the width dependence on the mass of tti at this 
particular value \jJ^ ax , for ir\ — > tttj and tti — > nrf respectively. The reason for this is to 
see the largest values one can obtain for these decays. It turns out that these values have 
a maximum on the order of 10 MeV. However, the value of should not differ much from 
/^tt and thus the above widths will be on the order of 1 MeV. 

In Tables 18.21 and 18.31 we presented the widths only for the modes including mesons 
with the lowest radial quantum number such as it and K. It allowed us to choose the free 
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Figure 8.13: Decay width for tt\ — ► nr] at fi v = [i™ ax as a function of the mass of tt\. 




Figure 8.14: Decay width for m — ► 7tt/ at fx v = \i™ ax as a function of the mass of n±. 

parameters \i equal to either /i^ or fix- This, however, cannot be justified for decays into 
mesons with higher values of the radial quantum number such as 77(1295) or p(1450). The 
dependence of the decay width on the value of /x^ for tt\ — > 717/(1295) with 5 = and m ex = 
2 GeV is shown in Fig. 18.151 This mode may be a significant one (the width is at most 
about 10 MeV), which agrees with the results obtained in jSHj. Unfortunately, there is not 
enough experimental data to constrain the free parameters /i for 77(1295) and other radially 
excited mesons. 

From these results it is clear that fully relativistic values are significantly different from 
nonrelativistic ones. There are two sources of this: the Wigner rotation which introduces 
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Figure 8.15: Decay width for tt\ — > 7r?7(1295) as a function of the parameter for a 7Ti 
with mass 2.0 GeV. 

relativistic coupling between spin and spatial degrees of freedom in the wave function, and 
different relations between energy, momentum and the invariant masses (in the phase space 
and most of all in the orbital wave functions). Both corrections actually appear to introduce 
corrections as large as 10% and thus should be included in phenomenological models. 

We also observe that in our relativistic constituent model the S + P selection rule, which 
was mentioned in Chapter 2, is obeyed. According to this rule the favored modes for the 
7Ti decay are nbi, vr/i^, V a i,2, and both channels KK\. Our results support this prediction 
and agree with other nonrelativistic models |51| I54j . 
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Chapter 9 



Final state interactions 



In the previous chapters we constructed meson states and studied kinematical relativistic 
effects at the quark and gluon level. In this chapter we will estimate the size of corrections 
to the it i decays originating from the meson exchange forces between mesons. Since the 
particle number is not conserved and momenta of particles are on the same order as their 
masses, this problem should be treated in a relativistic formalism. We will begin with the 
Lippmann-Schwinger equation applied to the irp and nbi states. Then we will describe 
the computational procedure of solving the resulting integral equations. At the end of this 
chapter the numerical results will be given. 



9.1 The Lippmann-Schwinger equation 




Figure 9.1: Final state interaction nbi <-> wp. 
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Figure 9.2: The lowest order FSI corrections to the tti decay from the interaction irbi <-> 71"/?. 

Other models |511l54ll5^] predict larger widths for the process %\ —* 7T/9 than our model. 
The numerical values were given in Chapter 2. It is possible that this width is increased 
by final state interactions between the outgoing mesons. The b\ created in the process 
7Ti — > 7r&i can subsequently decay into it and oj, and then the u) can absorb the other ir and 
create the p, as shown in Fig. 19.11 Of course the reverse irp — > irbi process can occur as well. 
An uj exchange may occur more than once and we must sum up all possible amplitudes; the 
lowest order diagrams are presented in Fig. 19.21 A bold horizontal solid line represents a 
hybrid and normal horizontal solid lines refer to mesons. A vertical dashed line corresponds 
to a single ui meson exchange. In order to describe a total contribution of the final state 
interactions to the decay widths of tti, we need to solve the Lippmann-Schwinger equation 

T = V + VGT, (9.1) 

which is equivalent to summing over all diagrams, including those with intermediate hybrid 
states. 

Let the iri state be denoted by index a, and the states irp and irbi by Roman letters. 
For simplicity we will assume the iri spin wave function has only the S = qqg component; 
that gives the largest widths for both ixb\ and up modes. Thus our calculations should give 
the upper limit for FSI corrections. Let us introduce the matrix elements for the potential 
V: 

V ai = (a\V\j), V tj = (i\V\j), (9.2) 

and similarly for T. The elements V a i are just the amplitudes of the corresponding decays 
of 7Ti, whereas are related to the final state interaction potential. In our state-space the 



Figure 9.3: The Lippmann-Schwinger equation shown diagrammatically. 
Lippmann-Schwinger equation can be written as 




Vji + VjkGkTki + Vj a G a T a i, 



(9.3) 



where 



Gi(E) = [E-H (i) + ie] 



-l 



(9.4) 



Both equations (|9.3|1 are represented diagrammatically in Fig. 19.31 A bold vertical solid line 
corresponds to the total amplitude of the interaction between two dimeson states, whereas 
a normal vertical solid line represents the sum over all intermediate states. 

Hereinafter, if an index appears twice or more, then summation over this index is im- 
plicitly assumed. The Hamiltonian of the state \i) is Hq(i), e — ► + , and E is the energy 
which will equal the mass of the tt\. Calculating TJj from the second equation in (|9.3|) and 
its substitution to the first one leads to 



Now we introduce the T matrix acting only between states \irp) and |7r&i), denoted by t 
and defined by 



Diagrammatically this represents the sum of all diagrams without hybrid intermediate 



Toi = V a i + V a jGj[l — VG\- k (Vki + Vk a G a T a i). 



(9.5) 



t = V + VGt, % = [l-FG]r fe %, 



(9.6) 
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Figure 9.4: The Lippmann-Schwinger equation in a meson subspace. 
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Figure 9.5: The self-energy of a hybrid. 



states, as shown in Fig. 19.41 Therefore 



T a i — V a i + V a jGjtji + V a jGj[l — VG]jfcVk a G a T a , 

= Vai + V a jGjtji + V a jGj[5jk + (tG)jk] Vk a G a T a {, 



(9.7) 



where we used [1 — VG] 1 = 1 + tG in the last term. 

Using the identity [1 — GV]™ 1 = 1 + Gt on the first two terms on RHS leads to 



T Qi = (V[l - GV]~ l ) m + [V(l + Gt)GV] aa G a T a 



(9.8) 



Solving for T a { gives 

T ai = [1 - (V(l + Gt)GV) aa G a ]-\V[l - GV}- 1 ^ = G" 1 ^ 1 - E a ]-\V[1 - GV}- 1 )^, 

(9.9) 

where S Q = (V(l + Gt)GV) aa is the self-energy of the tt\, shown in Fig. 19.51 Both G a 
and E„ are numbers, and thus 



T ai = -^(V + VGt) ai , 

hi — m ex — Li a 



(9.10) 



no 



Figure 9.6: Renormalized FSI correction to the tt\ decay amplitude. 

where we again used [1 — GV] -1 = 1 + Gt. Renormalization of the theory is related to 
cutting off the self-energy term. Therefore we obtain finally 

T(xi — Vai ^ajGjtji, (9.11) 

with tji defined by ()9.6|) . This is equivalent to excluding an exotic intermediate state from 
the FSI corrections. Diagrammatic representation of Eq. ()9.11|) is given in Fig. 19.61 

Eq. (|9.6I) . written in a full notation in the rest frame of the 7Ti, is 

t(p, q, A, A') = V(p, q, A, A') + / , , 3A jr. Ha V &> k ' A ' A " )G(fc)t(k, q, A" , A'), 

(9.12) 

where 



Ui(k) = E(rrii,k), H (k) = ui(k) +u 2 {k), 

G(k) = [E-H (k)+ie}-\ (9.13) 

and rrii (i = 1, 2) are the meson masses in an intermediate two-meson state which is related 
to G. In the above p and q are the relativistic relative momenta between two mesons, and 
the dependence on the center-of-mass momentum has been already separated. Spins A refer 
to either b\ or p. 

We introduce the partial wave potentials 

V L uM= E f dn p dn n (L,M;l,X\J,j)(L',M';l,X'\J,j) 

M,M',\Xj 



r\r\ 



xV(p, q ,\,\')Y LM (p)YZ, M ,(q), 



(9.14) 



where dQ^ is the element of the solid angle in the direction of the vector k and k = |k|. 
Similarly we define ti,y(p,q). Substitution of VLL'(p,q) into (|9.12|) gives 

C k 2 dk 

t L L>{p,q) = V LU (p,q)+J2 / TTTTTa TU 77xVLL»(p,k)G L »(k)t L » L ,(k,q), (9.15) 

^fj (2vrf 4w I /> i i(£;)u; L " i 2(£:) 

where 

w L »,i(fe) = E(m L » ti , k), Gl" (k) = [E — uj v%1 {k) - oj L " <2 (k) + ie]" 1 . (9.16) 

In our state-space we can have L = 0, 2 (the relative angular momentum between tt and 
&i) or L = 1 (between tt and p). For a 7Ti we also have J = 1. Thus in Eq. (|9.15j) we must 
substitute: 

mL,i = m n , mo i 2 = m 2) 2 = m i)1 , m\,i = m p . (9.17) 

From the conservation of parity the only nonzero components of the final state interaction 
potential are Voi, V\q, V\i and V21. Moreover, from the CP conservation for strong inter- 
actions we have Voi = V^q and Via = ^a* • The integral equation (|9.15|) cannot be solved in 
general analytically and one need to replace it with a set of matrix equations. The details 
will be given in the next section. When this is done and t^y are found, we may go back to 
Eq. (|9. llf) which becomes 

r k 2 dk 

tun = -up) + e / (ar) .(2j + i)^(^ w ^w G '^ t "'*- n (9 - 18) 

where cll{P) are the partial amplitudes defined in ()8.13j) and l|8.21|) . In order to obtain the 
corrected widths, these amplitudes must be replaced by the corrected amplitudes cll(P) in 
Eq. <H3EJ>- 

Finally, we proceed to the form of the final state interaction potential. The Lagrangian 
of the piriv vertex is 

Lpw = gpn^^d^dxpi, (9.19) 



and for the biiruj vertex is 

£&ittw = 9b 1 iTLud^\ % u^ (9.20) 

where g p7Til j and gbmoj are the hadrodynamical coupling constants and the index i corresponds 
to isospin. The first one lies between 0.01 MeV" 1 and 0.02 MeV -1 and we will take a value 
0.014 MeV -1 |79j . whereas the second one can be obtained by calculating the width of the 
decay b\ — > txlo and comparing with its experimental value 142 MeV [ZE]- Consequently, 
the amplitude of b\ — > ixuj is equal to 

A(P,A fel ,A w ) = - 5f)1 ™e*(A fel K*(A w ,P), (9.21) 

where P is given by 

E(m n ,P)+E(m u! ,P)=m bl . (9.22) 
The width for the angular momentum L between tt and to is 



P\Af 
8vr(2L + l)m 2 bi 



T L = _, nr \\, 2 , (9-23) 



where \ A\ 2 is summed over A^ and averaged with respect to A& 1 . With the help of 

£ e"(A, p)e"*(A, p) = -gT + (9.24) 

we have 

l^| 2 =^(l + £^-), (9.25) 

and from T s - wave = 142 MeV- 0.92 = 131 MeV (the factor 0.92 comes from the partial wave 
D/S ratio) we get g^-nuj = 3650 MeV. Thus gFSi = gpnu • gtmuj is around 50. 

The final state interaction potential can be obtained from Lagrangians (|9.19|) and (|9,20|) . 
written in the momentum space and dressed with the instantaneous u> propagator, 

V(p,q_,X bl ,X p ) = gFSie Pj varP^q v - ( To - o"e CT (Afe 1 ,-p)e T *(A p ,q), (9.26) 

(p - q) + m i 

where p is the momentum of the ir in the |7r&i) state and q is the momentum of the p. For 



large values of p and q this potential grows with no limit which is a consequence of treating 
mesons as elementary particles. Therefore we must regulate this potential with an extra 
factor that tends to zero for large momenta. We choose an exponential function, 



of the quantity which is invariant under translations and rotations. Here A is a scale 
parameter of order 1 GeV. This scale estimates a limit of treating mesons as elementary 
particles interacting via Lagrangians ()9.19j) and ()9.2Uj) . 

9.2 Computational procedure 

The Lippmann-Schwinger integral equation ()9.15j) corresponds to outgoing wave boundary 
conditions. This means that the singularity of the term G(k) is handled by giving the 
energy E a small positive imaginary part ie. An integral of this form may be solved using 
the Cauchy principal-value prescription, 



e 



|P*|/A 



(9.27) 




(9.28) 



where the principal value is defined by 




(9.29) 



with ko being the zero of the real function f(k). From the formula 




(9.30) 



one obtains 




(9.31) 



that can be generalized to the so-called Hilbert transform of a function /: 



P 



f(k)dk 



[f(k) - f(k )]dk 



U2 _ U2 



(9.32) 



The integral in Eq. (|9.15|) has a slightly different singular denominator term [E — Ho{k)\ , 
and in this case 



f(k)dk 



o E-H {k)+ie P J E - H (k) 

f°° f{k)dk k 2 -k 2 Q . fn .LOi(k )uj 2 (ko) 



f(k)dk . ( dHoW y 1 



k=ko 



Jo k 2 - k 2 E - H (k) 
oo dk rf/n/1.2 1.2 



^-^2 
00 d& 



f(k)(k' 2 - k'o) _ A; 2 fc o 1 

£ - flb(fc) fc ^ fc o J5 - Ho{k) 

f(k)(k 2 -k 2 ) 2f{ko)u 1 {ko)u 2 {ko) 
E - H (k) E 



^l(fcp)^2(feo) 

koE 



. u;i(fc )q; 2 (fco) 

«tt/(A;o) — . 9.33 

koE 



Making a transition 



/(*) 



k 2 V(p,k)t(k,q) 



(9.34) 



4(27r) 3 a;i(fc)a;2(fe)' 

and including the summation over partial waves, leads to a desired equation for t^yip^q) 
that no longer has a singularity: 



f°° dk 
tLV (p, q) = V LL , (p, q) + £ ^ 4(27r) 3 (fe2 _ A; 2 i „ ) 

H ^r- Vll" (p, A; ,l" )*l"L' (&o,l" , g) 



L" 



k ( k ~ k 0,L") V LL"(p, k)t L » L ,(k,q) 

[E - H ' L >>(k)]wL'\i(k)wL'\2(k) 
k o,n> V LL » (p, k , L " )twv (k),L» 



4(2vr) 3 £ 



-, (9.35) 



where 



H ,L(k) = u Li i(k) + u Lj2 (k) 



(9.36) 



and the quantities &o,L are defined by 



Ho, L (ko, L ) = E. 



(9.37) 



Because in the decay tv% — * nbi the D-wave amplitude is negligible as compared to that 
of the S-wave, we may reduce our angular momentum space to L = 0, 1. Therefore we need 



-i no 



only the expression for Vol ip, q) ■ From the definition (|9.14l) and with the help of 1)5.31(1 we 
obtain 

£ J dn p dn (l V(p,q,X,X')e ljk e i *(X)e^X')q k /q, (9.*S) 



Voi(p, q) 



iV3 



47T 



A. A' 



where V(p, q, A, A') was given in (|9.26[) . The S-matrix is obtained from the t matrix via 



x .y/k ,Lk ,L' , . 

JLL' = OLV - I — lQ n 2E t LL'\^0,L, ko,L')> 



(9.39) 



and for a 2 x 2 unitary matrix may be parametrized by scattering phase shifts 5q and 8\, 



S 



LV 



rje 



2i5\ 



(9.40) 



with < rj < 1. 

Eq. (|9,35|) may be solved, as we already mentioned, by converting the integration over k 
into the sum over N integration points k n , n = 1, 2, ... N (determined by Gaussian quadra- 
tures) with weights w n |8Uj : 



k l{kl - fcg L »)V LU '(p,k n )t L » L >(k n ,q) 



E 



Vll" ip, fco.L" )twv (ko,L" , q) 



L" 



[E - #0,L"(^n)]^L",l(&nVL» 5 2(fc 
fep,L" VlL" (gi fcp.L" )*L"L' (fep,L" , g) 



4(2vr) 3 £ 



(9.41) 



By taking p, q equal to either k n or fcp^ we can rewrite this equation as the (2N+2) x (2N+2) 
matrix equation 

tmn = Vmn + Rmoton, m,n,0= 1..2N + 2, (9.42) 

where 



V mn = to = 1..N + 1, n = l.JV + 1; 

Knn = V Q1 (l m ,l n ) m = 1..7V + l, n = A^ + 2..2AT + 2; 

Knn = %(im,in) to = A^ + 2..2AT + 2, n = l.JV + 1; 

V mn = to, = AT + 2..2N + 2, re = A r + 2..2AT + 2 



(9.43) 



and 



l n = k n n = l.JV; 
In+1 = ^0,0; 

l„ = k n - N - 1 n = N + 2..2N + 1; 
I2N+2 = ko,i- 



(9.44) 



The matrix R is defined by 



(9.45) 



with no summation over n, where 



D 



k±w n 



n 4(2n)*uJo,i(ln)u>o,2(l n )lE-Ho,o(ln)} 



D 



ko,0 



N + 1 — 4(27r) a E 

D 



TV 2k 0>o w n 
3=n kl-k\ 



2 

n "'0,0 



ktW n 



11 4(27r) 3 CJl,l(inVl,2an)[S-^0,l(in)] 



£>2AH 



4(27r) 3 _B 



2^n=l k 2 -k 2 



n = l.JV; 



n = N + 2..2N + 1; 



(9.46) 



Having inverted i? mn and solved for t mn , we subsitute it into (|9.18|) that in 2N + 2- 
dimensional space has the form 

N+l 2N+2 

«o(^0,o) = °o(^0,o) + X! Ai°oOn)*n,iV+l + X/ Ai°lOn)*n,iV+lj 
n=l n=JV+2 
JV+1 27V+2 

5i(A;o,i) = °i(%i) + 2J D nao{ln)t n ,2N+2 + ^2 D nai{ln)t n ,2N+2- (9.47) 
n=l n=Af+2 

The quantities &q o and fc^i are according to (|9.37j) the relative momenta of trbi and irp 
systems, respectively. Thus Eq. (j9.47j) gives the FSI-corrected amplitudes of tt\ — > 7r6i and 
7Ti — > 7T/9 as functions of the energy E that is equal to m ex . The S-matrix in terms of t mn is 



_ - y H(N+l)lj(N+l) 



i(N+l),j(N+l), 



(9.48) 



where i,j = 1,2 refer to channels 7r6i and 7rp. 



-i nr 



9.3 Numerical results 



The partial wave potential Q9.38JI obtained from the nbi — Trp final state interaction potential 
(j9.2fij) is not separable, i.e., cannot be represented as a product of functions of only one 
variable, 

V(p,q)^f(p)g(q). (9.49) 

This corresponds to non-locality of this potential in the position space and results from a 
finite structure of mesons. If they were elementary and therefore V(p, q) was a separable 
potential, then Eq. (|9.35|) could be solved analytically. In this section we will assume that 
this is the case and use it to test the accuracy of numeric results from the preceding section. 
Consider the following symmetric and off-diagonal potential, 



V L L'(p,q) -- 

We seek a solution of (j9,35|) in the form 

tLL'(p,q) = 



1 

1 



g{p)g{q)- 



( 



AoO Aoi 

V Aoi An 

Substitution of the above matrices into (|9.35|) leads to 



g(p)g(q)- 





Aoo 


Aoi 














1- 




:)•( 


V 


Aoi 


An j 









where the quantities Iq and I\ are given by the integrals: 

k 2 dkg 2 {k) 



h 



o 4(2vr) 3 u;o,i(£;K2 (*:)[£ " #o,o(fc) + 
k 2 dkg 2 (k) 

4(27r)3a;i i i(A;)wi,2(fc)[£; - H 0tl (k) + ie] 



(9.50) 



(9.51) 



(9.52) 



(9.53) 
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Figure 9.7: Phase shift Si for a toy interaction Trbi — irp as a function of the number of grid 
points N. 

For the t matrix we obtain then 



tLL'iP, q) 



J i 1 \ g{p)g{q) 



(9.54) 



and the phase shifts can be computed from (|9.39f) and 1)9.40(1 . The integrals Iq and I\ 
may be easily calculated using the principal-value prescription described in the preceding 
section. 

In Fig. 19.71 we show the numerical values of the phase shift 5\ in degrees as a function 
of the number of grid points N, and compared to the corresponding analytical value. The 
toy potential was chosen such that 



g(p)g(q) = -v e~^ + ^ 2 



(9.55) 



with Vq =1000 (the FSI potential is dimensionless) and /x = 1 GeV, whereas the energy 
E = m ex = 1600 MeV. We see that the numeric results for this phase shift converge pretty 
fast to a value that is in a good agreement with the analytical one. Accurate and stable 
numbers are obtained already for iV being approximately 10. 

Now we move to the original final state interaction and replace our toy potential with 
the real potential (|9.38j) obtained from Eq. (|9.26|) . Fig. 19.81 and Fig. 19.91 show an increasing 



stability of phase shifts <5o and Si (in degrees) for gFSi = 200 as the number of grid points 
N grows. The default values in these plots are E = 1600 MeV and A = 2 GeV. We see that, 
for an unseparable potential, a larger TV" is needed to obtain stable results. For N = 20 the 
discrepancy is still about 1°. The unitarity condition is satisfied here with accuracy better 
than 0.1%. 
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Figure 9.8: Phase shift So for the interaction irbi — irp as a function of the number of grid 
points N, for energy E = 1.6 GeV and gFSi = 200. 
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Figure 9.9: Phase shift Si for the interaction irbi — np as a function of the number of grid 
points N, for energy E = 1.6 GeV and gFSi = 200. 

Let us assume that the coupling constant gFSi is the variable and we introduce the 
potential strength I = gFSi/ 9fsii where gpsi = 50. In Fig. 19.101 and Fig. 19. Ill we present 



-1 nr> 



the phase shifts as functions of I. Their dependence on the energy E for I = 1 is shown in 

Fig. EH 




Figure 9.10: Phase shift 8q for the interaction irbi — up as a function of the potential 
strength /, for energy E = 1.6 GeV. 
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Figure 9.11: Phase shift 5± for the interaction irbi — np as a function of the potential 
strength /, for energy E = 1.6 GeV. 

Within an experimental range of the hybrid mass m ex = E and the hadrodynamical 
coupling constant g the phase shifts for both channels L = (7r&i) and L = 1 (np) are 
rather small. Therefore one expects that final state interaction would not change much the 
widths of 7Ti — > nbi and tt\ — > 7T/9. 

In Tables I97T1 and 19^21 we compare the original widths (in MeV) with the FSI-corrected 
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Figure 9.12: Phase shifts 5q and 5\ for the interaction irb\ — irp as functions of the energy 
E, for 7 = 1. 

ones for various values of m ex and A, for 1 = 2. The strength was doubled in order to 
recompense damping caused by the regulating factor in the FSI potential. For A other 
than 2 GeV, the coupling constant gFSi was renormalized such that the value of Voi(po ; Qo) 
remained equal to that for A = 2 GeV. Here po and qo are the physical values of the relative 
momenta between two mesons in the irb\ and irp channels, respectively. In Fig. 19.131 and 
Fig. 19.141 we show the width dependence on / for m ex = 1.6 GeV and A = 2 GeV. The 
behaviour of 5\ indicates the existence of the irp resonances for the free parameters given 
(a resonance occurs when the phase increases over n ■ 90°, where n = 1,2,3...). 

The numerical results confirm our predictions. Only large FSI potentials may change the 
width values for both modes by more than few MeV. However, the up channel is subjected 
to a rather large relative correction. It may be caused by a big value of the original width 
for the irb\ mode. In any event, the width for the ir\ — > irp is always on the order of a few 
MeV which agrees with [HI] or [HSJ and does not, for example, with |81j . 

There is one constraint for the possible values of the cut-off parameter A and the coupling 
constant gFSi-, given by the decay of the 7T2(1670). This meson almost does not decay 
into irbi but significantly does into irp. If the corrections from the final state interaction 
irb\ *-* irp were large, the irb\ mode width could be increased by the irp channel, resulting 
in a disagreement with experiment. Thus we arrive at the conclusion that FSI cannot 
introduce corrections that are too large. This agrees with our predictions for real I. 
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Table 9.1: Original and FSI-corrected widths in MeV for the irbi mode of the 7Ti decay, for 
various values of m ex and A. 
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Table 9.2: Original and FSI-corrected widths in MeV for the irp mode of the tt\ decay, for 
various values of m ex and A. 




Figure 9.13: FSI-corrected width of n\ — > irbi as a function of the potential strength /, for 
a 7Ti with mass 1.6 GeV. 




Figure 9.14: FSI-corrected width of 7Ti — ► 7rp as a function of the potential strength I, for 
a 7Ti with mass 1.6 GeV. 



Chapter 10 



Conclusions 



The quantitative description of the dynamics of gluons is complicated by the strong char- 
acter of their interactions. Hadrons with excited gluonic degrees of freedom seem to be the 
only way to obtain information about the nature of low-energy gluons. Therefore they are 
of a great importance for our understanding of the quark-gluon interaction. In this work 
we focused on exotic mesons because they give insight into the dynamics of normal mesons. 

The current candidates for the lightest exotic meson are 7Ti(1400) and 7Ti(1600). Most 
of the experimental resonances have large widths compared to normal mesons that decay 
strongly. Moreover, the observed decay channels should not be dominant according to the 
theory. Therefore, the reported signals do not completely agree with theoretical expecta- 
tions and may originate from different phenomena such as rescattering. In order to resolve 
this problem it will be necessary to observe and measure exotic states with other quantum 
numbers. Searches for these states, as well as further exploration of the tt\ decays, are 
planned at Jefferson Lab (GlueX). 

The dynamics of exotic mesons at the scale of ~ 1 — 2 GeV should be described by a 
relativistic theory. In this work we studied the size of relativistic effects in the decays of 
the 7Ti, and discussed a new picture of meson decays in the Coulomb gauge. Our considera- 
tions led to two important conclusions. First, numeric results showed significant relativistic 
corrections arising from the spin-orbit correlations introduced by Wigner rotation. The 
widths calculated using fully relativistic formulae are in general larger than the correspond- 
ing values calculated with no Wigner rotation (by a factor of order 10%), but smaller than 



completely nonrelativistic values. Some decays that are suppressed in the nonrelativistic 
limit (for example tt\ — > up if the orbital wave functions of ir and p are the same) acquire 
nonzero amplitudes in the relativistic case. 

Second, the tti prefers to decay into two mesons, one of which has no orbital angular 
momentum and the other has L = 1 (the S + P selection rule). Thus this selection rule, 
found in other models seems to be quite general. Some decays (71-77, P 1 ^-, KK*) are suppressed 
by symmetries in the orbital wave functions, and the assumption that the parameters p for 
mesons with the same radial quantum numbers should be almost equal. We also noticed 
that the rates for the higher partial waves in decays of normal mesons where two waves are 
possible are larger in the 3 Po then the 3 Si model. 

However, there is one problem which needs to be explored further. There are several 
components of the qqg normal and exotic meson wave functions, and in order to give accurate 
numerical predictions for the widths, one needs to know the relative amounts of these 
components. Since we were interested in relativistic effects and not in the absolute width 
values, we calculated the amplitudes for each component separately. 

We also studied the interaction between irbi and irp final states. For the physical range 
of the parameters used in our analysis we observed rather small absolute corrections arising 
from FSI. The decay width for the process 7Ti — ► irbi decreased only by a factor of order 1 
MeV, whereas for the it\ — > irp increased by a similar amount. Therefore, the latter process 
seems to be suppressed anyway. 

The experiments have reported exotic resonances that are light and rather broad. This 
work showed that relativistic calculations give widths smaller than nonrelativistic ones, thus 
the discrepancy between theory and experiment becomes larger. For the dominant 7r6i mode 
of the 7Ti(1600) decay the maximum width is ~ 200 MeV as opposed to the experimental 
~ 300 MeV. Furthermore, the width in the irp mode is on the order of a few MeV, and the 
modes irrj and tttj' are negligible. 

The experimental data for the tt\ exotic meson are not well understood yet. Further 
experiments are needed to clarify the nature of the reported signals and explain the dynamics 
of exotic mesons. This work showed that, in order to compare these data with the theory, 
one should use models which are relativistic. 
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